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HOT  FILM  VELOCITY  MEASUREMENTS  DOWNSTREAM 
OF  A  SWEPT  BACKWARD-FACING  STEP 

Chapter  1 
INTRODUCTION 

Overview 


Over  the  past  25  years  many  investigators  have  been  interested  in  the 
flow  over  an  unswept  backward-facing  step.  Through  their  experimental  and 
computational  efforts,  a  modest  attempt  has  been  made  at  understanding 
the  resulting  flowfield.  Because  of  the  large  niimber  of  parameters  required 
to  accurately  define  the  flowfield  and  reattachment  point  downstream  of  a 
backward-facing  step,  the  work  has  progressed  slowly.  Even  today,  all  the 
parameters  still  have  not  been  identified.  Recently,  some  interest  has 
developed  in  the  more  general  flowfield  over  a  swept  backward-f<icing  step,  of 
which  the  unswept  step  is  a  special  case.  Very  little  experimental  data  and 
theoretical  analysis  exist  for  the  swept  backward-facing  step. 

The  purpose  of  this  study  is  to  provide  preliminary  measurements  of 
the  mean  velocity^  profiles  at  various  survey  locations  up-  zuid  downstream  of 
a  backward-facing  step  with  sweep  angles  of  0,  15,  and  30  degrees.  It  is  felt 
that  the  results  from  this  study  can  be  applied,  at  least  indirectly,  to  mziny 
of  the  problems  involving  this  type  of  flow.  Severed  examples  of  where  the 
flow  over  a  swept  backward-facing  step  could  be  applied  eire:  the  flow  around 
buildings,  cliffs  and  moimtains;  the  flow  through  diffusers  and  combustion 
chambers;  and  the  flow  over  seams  of  an  airplane  wing  or  over  the  rotating 
band  of  a  nutating  artillery  projectile.  One  very  important  use  for  the 
results  from  a  study  of  the  flow  over  a  swept  backward-facing  step  is  in  the 
development  and  verification  of  computational  turbulence  models,  both 
two-  and  three-dimensional. 

In  order  to  understand  the  flow  over  a  swept  backward-facing  step 
several  aspects  of  the  flowfield  mxist  be  defined,  such  as:  the  initial  conditions 
at  the  point  of  flow  separation,  the  flow  in  the  recirculation  region,  the  flow 
reattachment  region  and  the  recovery  region  downstream  of  reattachment. 
Eaton  and  Johnston  (1981),  in  their  review  of  the  research  to  date  of  the 
flow  over  unswept  backward-facing  steps,  listed  several  reasons  why  this 
geometry  is  the  best  choice  for  investigating  turbulent  flowfields.  First, 
although  the  flowfield  itself  is  extremely  complex,  the  flow  downstream  of  a 
backward-facing  step  is  one  of  the  simplest  two-dimensional  reattaching 
flows  known.  The  flow  sepsirates  along  a  straight  line  represented  by  the 
edge  of  the  step.  At  flow  separation,  the  streamlines  are  nearly  peirallel  to 
the  wall.  Finally,  the  flow  within  the  recirculation  region  is  highly  turbulent 
and  the  boimdary-layer  is  separated.  The  same  reasons  can  also  be  applied 
to  the  three-dimensional  flow  over  a  swept  backward-facing  step. 

Research  Objectives 

In  this  study,  three  0.5  inch  high,  backward-facing  step  configurations 
with  sweep  mgles  of  0,  15  and  30  degrees,  respectively,  were  tested  in  a 
subsonic  wind  tuimel.  A  step  height  of  0.5  inches  was  selected  in  order  to 
allow  comparisons  to  be  made  with  Selby’s  (1982)  work.  Selby  performed 
numerous  flow  observations  over  a  0.5  inch  swept  backward-facing  step.  Oil 
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flow  observations  were  performed  in  order  to  determine  if  the  flowfields  were 
similar  to  those  obtained  by  other  investigators.  A  hot  film  anemometer 
system  was  used  to  measure  the  components  of  the  total  velocity  parallel 
and  perpendicular  to  the  step  face.  The  intent  of  this  study  was  to  provide 
preliminary  data  on  how  the  mean  velocity  was  affected  by  step  sweep  angle, 
so  that  a  more  detailed  survey  of  the  flowfield  could  be  accomplished  during 
future  testing.  Velocity  profile  surveys  were  conducted  at  locations  up-  and 
downstream  of  the  step.  Probe  positioning  was  a  manual  operation  with 
data  acquisition,  storage  and  reduction  being  computer  controlled.  From  the 
hot  film  measurements,  the  total  velocity  and  flow  angle  were  determined. 
Knowing  this  information,  the  velocity  components  parallel  and 
perpendicular  to  the  step  face,  as  well  as  an  estimate  of  the  bound2Lry-layer 
thickness  at  each  survey  location  were  calculated.  The  results  from  this 
study  were  also  applied  to  two  computational  turbulence  models.  The  first 
method  examined  was  Coles’  law  of  the  wall  and  wake  for  approximating 
attached  turbulent  boimdary-layers.  The  second  method  was  a  modified 
Baldwin-Lomax  algebraic  model  for  attached  and  sepau'ated  turbulent  flows 
over  a  backward-facing  step. 

Description  of  Flowfield 

The  flow  over  a  backward-facing  swept  or  unswept  step  represents  a 
relatively  simple  reattaching  flow  but  the  flowfield  itself  is  extremely 
complex.  For  the  unswept  step,  the  flow  is  basically  perpendicular  to  the 
face  of  the  step.  It  is  believed  that  the  flow  over  a  swept  backward-fcw:ing 
step  is  similar  to  the  flow  over  a  yawed  cylinder  as  described  by  Schlichting 
(1979).  By  using  this  comparison  the  flowfield  can  be  greatly  simplified. 

The  flowfield  is  no  longer  truly  S-dimensioned.  The  flow  over  an  infinite 
swept  step  is  then  composed  of  two  components,  one  parallel  and  the  other 
perpendicular  to  the  face  of  the  step.  The  perpendicular  component  is 
independent  (uncoupled)  from  the  other,  with  the  perpendicuW  flow  being 
analogous  to  the  flow  over  an  unswept  step  and  the  parallel  flow  being  that 
part  of  the  flow  caused  by  the  sweep  angle  of  the  step. 


A  diagram  of  the  different  regions  of  the  flowfield  over  a  swept  or 
unswept  backward-facing  step  is  illiistrated  in  Figme  1.1.  Above  the  surface 
of  the  step  (y  >  5),  within  the  freestreaun,  the  flowfield  is  inviscid.  Ordinary 
boundary-layer  growth  is  experienced  upstream  of  the  step,  with  this  part  of 
the  flow  having  little  effect  on  the  flow  downstream  of  the  step.  The 
upstream  flow  is  used  to  determine  the  fluid  properties  at  separation.  The 
edge  of  the  shear  layer  separates  the  inviscid  flow  of  the  freestream  from  the 
boimdctry-layer  flow  where  viscous  and  turbulent  stresses  dominate.  Just 
downstream  of  the  step  is  the  recirculation  region.  The  dividing  streamline 
separates  the  recirculation  region  from  the  rest  of  the  shear  layer.  The 
dividing  streamline  begins  at  the  edge  of  the  step  and  continues  downstream 
nearly  parallel  to  the  surface  of  the  wall.  Then,  as  the  reattachment  region 
is  approciched  the  dividing  streamline  curves  sheuply  down  to  the  surface. 
The  recirculation  region  is  below  the  dividing  streamline.  Within  the 
recirculation  region  there  are  several  areas  of  interest.  There  is  believed  to 
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be  a  comer  eddy  in  the  area  where  the  step  face  meets  the  lower  plate 
silrface.  The  major  area  of  interest  is  the  recirculation  region.  This  airea  is 
split  roughly  in  half  by  a  stagnation  line.  Above  this  line  the  fluid  in  the 
recirculation  region  is  in  the  direction  of  the  freestream  but  below  the  line 
the  flow  direction  is  reversed.  Within  the  recirculation  region  the  boundary- 
layers  are  separated  from  the  wall.  The  reattachment  region  is  located  at 
the  downstream  edge  of  the  recirculation  region.  Downstream  of  the 
reattachment  region  is  the  recovery  region  where  the  flow  is  realigning  itself 
in  the  direction  of  the  freestream.  Another  part  of  the  flowfleld  not 
indicated  in  Figure  1.1  is  the  crossflow  which  is  that  part  of  the  flowfleld 
parallel  to  the  step  face.  Note  that  with  parallel  flow  along  the  swept  step, 
the  perpendicular  flow  stagnation  line  is  not  necessarily  a  line  of  zero 
velocity. 

The  terminology  and  variables  used  to  describe  the  flow  over  a 
backward-facing  step  are  shown  in  Figure  1.2.  The  x  axis  is  along  the  wind 
tunnel  centerline  in  the  direction  of  the  free  stream.  The  y  axis  is 
perpendicular  to  the  wind  tunnel  floor  and  the  z  axis  is  in  the  direction  of 
the  tunnel  side  wall  with  the  positive  sense  to  the  left  when  looking 
downstream.  The  x'  and  z'  directions  are  perpendicular  and  parallel  to  the 
step  face,  respectively.  The  y'  axis  is  in  the  same  direction  as  the  y  axis. 
The  u  and  w  velocity  components  are  along  the  x'  and  z'  axes,  respectively, 
with  the  positive  sense  indicated.  The  step  angle  is  denoted  by  A.  The  total 
velocity,  U  and  flow  direction,  a  are  shown  with  their  positive  sense 
indicated. 


Literature  on  Swept  and  Unswept  Backward-Facing  Steps 

The  amount  of  literatme  available  on  the  flow  over  backward-facing 
steps  is  extensive.  Since  several  investigators  (Eaton  and  Johnston  (1981), 
Selby  (1982),  and  Hartman  (1988)],  have  provided  detjuled  historical  reviews 
of  work  performed  on  this  type  of  flowfleld,  one  will  not  be  provided  here. 
Instead,  a  brief  discussion  will  be  presented  of  those  references  used 
extensively  in  this  study  as  well  as  recent  published  works  on  the  subject.  A 
complete  list  of  cited  and  imcited  references  is  contained  in  the  Bibliography. 

Experimental  Technique  and  Procedure 

Several  texts  and  reports  describing  hot  wire/fllm  techniques,  data 
acquisition  and  reduction  were  used  to  establish  the  test  and  data  reduction 
procedure  used  in  this  study. 

Lomas’s  (1986)  text  on  hot  wire  anemometry  was  an  invaluable 
source  of  theoretical  and  practical  information  on  the  use  of  hot  wir^fllm  to 
measure  air  velocities  and  flow  directions.  This  text  describes  the  effects  on 
the  probe’s  output  as  a  function  of  the  probe’s  yaw  and  pitch  angle  to  the 
flow.  Also,  effects  of  humidity,  temperature  differences,  heat  conduction  to 
probe  supports,  measurements  close  to  a  wall  amd  many  other  experimental 
concerns  are  discussed. 
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Several  Thermal  Systems  Incorporated  (TSI)  technical  bulletins  were 
used  in  this  study  (Thermal  Systems  Incorporated  TB5,  TB16  and  TB18). 
Technical  Bulletins  16  and  18  describe  several  techniques  that  could  be  used 
to  compensate  hot  film  measurements  for  differences  in  temperature  between 
the  calibration  results  and  the  test  data.  Technical  Bulletin  5  provides 
general  information  on  the  use  of  hot  wire/fltm  anemometers  for  velocity 
measurements. 

The  data  reduction  procedure  used  to  determine  the  flow  velocity  and 
direction  from  a  single  component  hot  film  anemometer  is  described  by 
Olivari  (1978).  For  a  specific  data  point,  this  procedure  requires  that 
readings  be  taken  at  two  or  three  yaw  angles.  Then,  using  the  velocity  and 
yaw  c^ibration  results,  the  total  velocity  and  flow  direction  can  be 
determined.  Another  technique,  described  by  Al-Beirutty,  Arterberry  and 
Gessner  (1988),  uses  a  single  component,  slanted  hot  wire  to  measiire  the 
three  velocity  components  as  well  as  the  sbc  Reynolds  Stress  terms.  This 
procedure  could  prove  extremely  useful  in  future  velocity  and  turbulence 
measurements  in  the  flowfield  aroimd  swept  backward-facing  steps. 

Experimental  Results 

The  results  from  Selby’s  (1982)  phenomenological  study  of  the 
subsonic  flow  over  a  swept  backward-facing  step  proved  very  helpful  with 
the  experimental  setup  for  this  study.  Selby  performed  numerous  oil  drop, 
smoke  wire  and  surface-mounted  tufts  flow  visualization  tests,  which  were 
tiseful  in  establbhing  the  degree  of  similarity  between  the  flowfields  of  this 
study  and  his.  He  measured  the  mean  velocity  using  a  hot  wire  anemometer 
and  also  measured  the  static  pressure  along  the  surface  of  the  model.  From 
his  observations  and  measurements  he  was  able  to  determine  the  flew 
reattachment  location.  He  examined  what  effect  step  height,  sweep  angle, 
base  geometry  and  end  condition  had  on  the  flowfield  and  the  flow 
reattachment  location.  Two  important  findings  came  from  his  work.  First, 
for  sweep  angles  less  than  38  degrees,  the  w  velocity  component,  that  part  of 
the  flow  parallel  to  the  step,  could  be  uncoupled  from  the  two-dimensional 
flow  normal  to  the  step.  Second,  the  u  velocity  component,  the  flow  normal 
to  the  step,  behaved  like  the  flow  over  a  unswept  backward-facing  step. 

This  is  similar  to  the  flow  over  a  yawed  cylinder  described  by  Schlichting 
(1979). 


The  experimental  work  performed  by  Bradshaw  and  Wong  (1972) 
provided  data  for  comparison  with  the  mean  velocity  flow  over  the  unswept 
backward-facing  step  near  reattachment. 

A  good  source  of  experimental  data  and  an  excellent  historical  review 
of  research  on  the  flow  over  a  backward-facing  step  was  presented  by  Eaton 
and  Johnston  (1981).  Eaton  and  Johnston  disciissed  some  of  the  parameters 
that  affect  the  flowfield.  These  parameters  include:  the  initial  boundary-layer 
state,  initial  botmdary-layer  thickness,  freestream  turbulence,  pressure 
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gradients  and  the  ratio  of  step  height  to  channel  width  (aspect  ratio,  AR). 

A  recent  study  performed  by  Ellzey  and  Berbee  (1988)  investigated 
the  effects  of  Reynolds  number  and  aspect  ratio  on  the  flow  behind  a 
backward-facing  step.  They  foimd  that  there  is  little  effect  on  the  mean 
velocity  profiles  down  to  an  AR=4,  which  is  considerably  lower  than  the 
value  of  10  recommended  by  deBr(^erode  and  Bradshaw  (1972). 

Computational  and  Theoretical  Literature 

Compared  to  experimentaJ  studies  there  are  far  fewer  computational' 
and  theoretical  investigations  and  most  of  these  have  been  limited  to  the 
flow  over  the  unswept  backward-facing  step. 

A  recent  investigation  by  Awa,  Kline  and  Ferziger  (1988)  for  the 
flow  over  an  unswept  backward-facing  step  employed  a  zonal  approach  in 
determining  the  turbulent  flowfleld.  Their  method  used  a  second-order 
accurate,  differencing  scheme  applied  to  the  diflerent  zones  of  the  flow 
behind  the  backward-facing  step.  Their  technique  yielded  a  better 
prediction  of  the  turbulence  terms  than  could  be  obtained  from  the  standard 
k-e  method. 

A  finite  element  analysis  was  conducted  by  Thomas,  Morgan  and 
Taylor  (1980)  for  the  flow  over  a  unswept  backward-facing  step.  Their 
procedure  was  applied  to  both  laminar  and  turbulent  flows  and  achieved 
acceptable  agreement  with  experimental  results. 

Hu,  Wang,  Danberg  and  Seidel  (1984)  developed  a  finite  element 
method  for  calculating  an  incompressible,  laminar  flow  over  a  swept 
backward-facing  step.  Their  procedure  separated  the  flow  into  two  parts, 
the  first  being  along  the  step  face  and  the  second  perpendicular  to  the  step 
face.  The  two  flows  were  then  solved  independent  of  eaxh  other.  They 
indicated  that  their  results  looked  promising,  but  for  their  method  to  be 
useful  it  should  be  expanded  to  cover  tmbulent  flows. 

Work  performed  by  Hartman  (1988)  and  Hartman,  Seidel  and 
Danberg  (1989)  included  the  development  of  a  finite-difference  computer 
code  to  investigate  the  flow  behind  a  swept  backward-facing  step.  A  two- 
layer  eddy  viscosity  model  was  used  within  the  computer  code  which 
provided  for  the  presence  of  a  vertical  wall  and  flow  reattachment.  The 
computations  indicated  that  reattachment  length  decreased  with  sweep  angle 
and  step  height  and  these  findings  agreed  reasonably  well  with  experimental 
results. 


Several  investigations  which  proved  helpful  with  the  theoretical 
analysis  of  this  study  included  the  work  performed  by  Danberg  (1971).  In 
this  work  Danberg  describes  a  procedure  that  uses  a  least  squares  method  to 
determine  the  coefficients  of  an  equation  based  on  Coles’  law  of  the  wall  and 
wake.  This  procedure  can  then  be  used  to  approximate  the  velocity  profiles 
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of  an  attached  flowfleld. 

Danberg  and  Patel  (1989)  applied  the  Baldwin-Lomax  algebraic 
turbulence  model  to  the  flow  over  a  rotating  band  of  an  artillery  projectile. 
The  geometry  of  this  application  is  similar  to  a  forward  and  rearward  facing 
step.  A  detailed  discussion  of  the  Baldwin-Lomax  model  was  provided  along 
with  the  modifications  necessary  to  apply  it  to  the  flow  over  a  swept 
backward-facing  step. 

The  experimental,  computational  and  theoreticzJ  studies  cited  above 
are  not  intended  to  be  a  comprehensive  listing  of  all  work  on  the  subject  of 
the  flow  over  swept  and  unswept  backward-facing  steps.  Instead,  the  intent 
is  to  highlight  some  of  the  recent  developments  as  well  as  present  those 
investigations  which  pertained  directly  to  this  study. 
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Chapter  2 

EXPERIMENTAL  FACILITIES,  HARDWARE  AND  MODEL 

Overview 


This  study  was  conducted  at  the  U.S.  Army  Chemical  Research 
Development  and  Engineering  Center  (CRDEC),  Aberdeen  Proving  Ground, 
Maryland.  CRDEC  facilities  used  during  this  study  included  the 
Aer^ynamic  Research  Concepts  and  Assistance  Branch  (ARCA  Br)  and  the 
Experimental  Design  and  Fabrication  Branch  (EDF  Br|.  The  ARCA  Br 
provided  the  subsonic  wind  tunnel,  data  acquisition  ana  computer  systems; 
while  the  EDF  Br  built  the  wind  tiinnel  model. 

Description  of  Wind  Tiinnel 

The  continuous  flow,  open  circuit,  subsonic  wind  tunnel  used  for  this 
test  is  shown  in  Figure  2.1,  and  is  described  by  Miller  (1983).  Air  flow  in  the 
tunnel  is  created  by  a  125  horsepower  electric  motor  driving  a  squirrel  cage 
fan  located  downstream  of  the  test  section.  The  air  is  drawn  from  the 
outside  of  the  building  into  the  tunnel  settling  chamber  through  a  section  of 
honeycomb  and  two  fine  mesh  screens.  The  air  moves  from  the  settling 
chamber  through  a  6:1  area  contraction  into  the  test  section.  Dimensions  of 
the  tunnel  test  section  are  28  inches  high  by  40  inches  wide  by  88  inches 
long.  Downstream  of  the  test  section,  the  air  flows  through  a  diffuser,  passes 
the  fan  and  then  is  exhausted  through  the  roof  of  the  laboratory.  Tunnel 
velocit]'^  is  controlled  by  adjusting  a  series  of  vanes  in  front  of  the  fan.  The 
velocity  range  for  the  tunnel  is  approximately  10  fps  to  220  fps,  (0.009  to  0.2 
Mach  Number).  The  test  section  velocity  is  indicated  on  an  mclined 
manometer  attached  to  the  tunnel.  The  ttinnel  is  capable  of  continuoxis 
operation.  One  minor  problem  that  developed  during  testing  was  the 
i^uence  of  ambient  wind  gusts  on  the  test  section  velocity.  Wind  gusts  in 
the  direction  of  the  inlet  could  cause  the  test  section  velocity  to  fluctuate  as 
much  as  2  fps  during  a  test.  To  avoid  this  problem  the  outside  wind 
conditions  were  continually  monitored. 

Description  of  Traversing  Mechanism 

After  several  iterations,  a  Hnal  design  for  a  probe  traversing 
mechanism  was  settled  on.  Figure  2.2  shows  the  traversing  mechanism 
installed  in  the  roof  of  the  subsonic  wind  tunnel.  In  order  to  insure  accurate 
vertical  positioning  of  the  probe,  relative  to  the  surfeice  of  the  model,  a 
commercially  available  24  inch  height  gauge  was  8u:quired.  The  height  gauge 
scribe  was  replaced  with  a  specially  designed  adapter  which  allowed  the 
rotation  of  the  hot  film  probe  holder  about  the  longitudinal  axis  while,  not 
allowing  the  probe  to  move  vertically  towzo^ds  or  away  from  the  plate.  The 
adapter  incorporated  a  yaw  angle  indicator,  calibrated  in  one  degree 
increments  (0  to  360  degrees),  to  determine  probe  yaw  angle  or  flow 
direction.  A  thumb  screw  was  used  to  lock  the  probe  at  the  desired  yaw 
angle.  The  hot  film  probe  holder  rzui  through  a  thin  walled  steel  tube  which 
acted  as  a  wind  shield  and  both  were  secured  to  the  adapter.  The  wind 
shield  stopped  approximately  6  inches  above  the  probe  tip  so  that  it  would 
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Hind  Tunnel  Test  Facility 


Figure  2.2  Side  View  of  Traversing  Mechanism 
Installed  in  Wind  Tunnel 


not  adversely  affect  the  measurements.  A  slot,  milled  into  the  roof  and 
along  the  centerline  of  the  tunnel,  shown  in  Figure  2.3,  allowed  the 
traversing  mechanism  to  be  positioned  over  the  range  of  approximately  2.75 
inches  upstream  to  7,5  inches  downstream  of  the  step.  Two  guide  rods  were 
positioned  an  equal  distance  along  both  sides  to  keep  the  traversing 
mechamism  aligned  with  the  slot.  The  positioning  of  the  probe,  both 
vertically  and  longitudinally,  was  a  manual  operation.  For  future  testing,  it 
would  be  advantageous  to  automate  the  probe  positioning  operation  to 
reduce  errors,  fatigue  and  test  time. 

Description  of  Wind  Tunnel  Model 

The  wind  tunnel  model  was  fabricated  from  half  inch  plexiglass,  thus 
providing  a  half  inch  step  height.  Figures  2.4  through  2.6  give  the  overall 
dimensions  and  layout  for  the  backward-facing  step  wind  tiinnel  model.  A 
0.25  inch  aluminum  plate  was  used  to  support  the  base  of  the  model  and 
provide  for  a  solid  connection  to  the  floor  of  the  wind  tunnel.  Connection  of 
the  support  plate  to  the  tunnel  was  made  by  four  0.5  inch  threaded  rods 
located  near  the  support  plate  center  2nd  four  0.25  inch  tumbuckles  located 
at  the  comers  of  the  support  plate.  The  use  of  threaded  rods  and 
tumbuckles  allowed  for  setting  the  plate  at  an  angle-of-attack  to  the  flow  in 
the  wind  tunnel.  The  base  of  the  model  was  fabricated  from  0.5  inch 
plexiglass  and  was  designed  to  span  the  width  of  the  test  section;  thus,  the 
side  walls  of  the  test  section  acted  as  end  plates  for  the  model.  Li  order  to 
improve  the  flow  around  the  leading  edge,  the  base  plate  was  beveled  to  a  40 
degree  wedge.  The  overall  dimensions  of  the  base  plate  were  39.5  by  59.6 
inches.  The  final  components  of  the  model  were  the  three  interchangeable 
swept  backward-facing  steps.  The  first  step  was  unswept  or  had  a  0  degree 
sweep  angle,  the  other  two  had  a  15  and  30  degree  sweep  angle,  respectively. 
All  three  steps  had  a  centerline  distance  of  27  inches  measured  from  the 
leading  edge  to  the  step  edge.  It  is  important  to  note  that  the  upstrezun 
portion  of  each  step  was  constmcted  as  one  piece  so  that  there  would  be  no 
joints  between  the  leading  edge  and  the  step.  While  this  approach 
eliminated  edl  upstream  joints  it  presented  a  problem  in  that  there  were 
different  leading  edges  for  each  step.  Thxis,  the  upstrezim  flowfields  may  not 
have  been  identical  for  each  configuration. 

Initially,  the  leading  edge  of  the  three  step  configurations  were 
machined  to  continue  the  40  degree  wedge  from  the  base  plate  to  a  knife 
edge  along  the  upper  surfawre  of  the  step  plates.  As  will  be  explained  in  the 
next  chapter,  the  leading  edge  of  e2u:h  step  plate  was  radiused  in  order  to 
improve  the  flow  over  the  upper  surface.  All  model  hardware  were  secured 
together  with  screws  from  the  underside  so  that  the  upper  surface  remained 
as  smooth  as  possible.  Figtire  2.7  shows  the  0  degree  step  installed  in  the 
subsonic  wind  tunnel. 
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Figure  2.3  Top  View  of  Traversing  Mechanism 
Installed  in  Wind  Tunnel 
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Figure  2.4  Position  of  Model  Relative  to  Subsonic  Hind  Tunnel  Test  Section 


NOTE:  Small  circles  are  mounting  holes. 


All  diHcnslofis  ire  In  Inches. 


Figure  2.5  Model  Base  Hardware  (top)  0.25  Inch 
Aluminum  Plate,  (bottom)  0.5  Inch 
Plexiglass  Base 
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Two  physical  pzirameters  conunonly  used  to  describe  the  experimental 
apparatxis  and  setup  for  measuring  the  flow  over  a  backward-facing  step  are 
the  test  apparatus  aspect  ratio,  AR,  and  expansion  ratio.  Both  of  these 
parameters  affect  the  flowHeld  and  are  dependent  on  the  physical  dimensions 
of  the  wind  tunnel  and  model.  The  AR  of  the  test  apparatus  is  defined  as 
the  width  of  the  wind  tunnel  divided  by  the  step  height.  For  this  wind 
tunnel  facility  and  model,  the  AR  =  80.  If  the  AR  is  too  small  it  is  possible 
that  the  side  walls  of  the  test  section  could  adversely  effect  the  flow  over  the 
step,  de  Brederode  and  Bradshaw  (1972)  found  that  the  lower  limit  for 
ignoring  the  effects  due  to  AR  was  10,  but  Ellzey  and  Berbee  (1988)  foimd 
little  or  no  effect  down  to  values  as  low  as  4.  The  second  parzuneter,  the 
expansion  ratio,  is  defined  as  the  ratio,  relative  to  the  step,  of  the 
downstream  to  upstream  cross  sectional  areas.  The  backward  facing  step 
model  was  installed  in  the  wind  tunnel  so  that  the  upper  surface  of  the 
model  was  approximately  6.25  inches  above  the  test  section  floor.  The 
location  of  the  model  at  this  position  yielded  an  expansion  ratio  of  1.03.  The 
expansion  ratio  effects  the  flowfield  by  producing  a  streamwise  pressure 
gradient.  Kuehn  (1980),  showed  that  as  the  expansion  ratio  increased  so  did 
the  reattachment  lengtn.  Selby  (1982)  investigated  the  coupling  between  the 
wind  tunnel  model  and  the  test  section.  He  concluded  that  the  streamwise 
expansion  ratio  is  the  major  caxise  of  the  streaimwise  pressure  gradient  within 
the  separated  flow  region  of  swept  backward-facing  step  mounted  in  a 
typical  wind  tunnel  test  section.  A  compensating  roof  step  is  the  easiest  and 
best  method  of  eliminating  any  pressure  gradient  due  to  the  expansion  ratio. 
Since  the  expansion  ratio  for  this  study  was  very  close  to  unity,  and  any 
streamwise  pressure  gradient  was  assiuned  negligible,  compensating  roof 
steps  were  not  used. 

Description  of  Hot  Film  Anemometer 

The  velocity  profile  measmements  were  made  using  a  Thermo- 
Systems  Inc.  (TSI)  constant  temperature  hot  wire/film  anemometer  system. 
A  photograph  of  the  TSI  system  is  shown  in  Figure  2.8.  The  system 
consisted  of  a  1051-6  power  supply  and  analog  meter,  three  1050  Constant 
Temperature  Anemometers  (CTA),  and  three  1052  Linearizers.  For  these 
tests,  the  1052  Linearizers  were  not  used.  Typical  CTA  settings  were  as 
follows:  Bridge  Selection,  setting  1  under  5:1;  Probe  Resistance  Decades,  set 
to  manufacture  recommended  vdue;  High  Frequency  Filter,  set  to  100  kHz; 
ztnd  Reference  Voltage,  set  to  5  volts,  which  is  recommended  for  hot  film 
probes.  A  36  inch  probe  holder  provided  the  needed  capability  of  traversing 
from  the  ceiling  to  the  floor  of  the  wind  txmnel  test  section.  A  TSI  1210-20 
Hot  Film  single  component  probe  was  ;ised  to  measure  the  velocity  profiles 
throughout  the  experiments.  The  1210-20  Hot  Film  probe  has  a  sensor 
diameter  of  0.002  inches  zuid  a  length  of  approximately  0.040  inches.  These 
probes  also  have  a  recommended  minimum  and  maximum  velocity  range  of 
0.5  to  1000  fps,  respectively. 


18 


19 


2.8  Thermal  System  Incorporated  (TSI)  Hot  Wlre/FIlm 
Anemometer  System 


The  capability  to  measure  flow  angles  in  at  least  one  plane  and 
availability  of  hot  film  probes  were  the  main  reasons  for  the  selection  of  a 
hot  film  anemometer  to  measure  the  velocity  profiles  in  the  flowfield  around 
a  swepted  backward-facing  step.  Although  other  instrumentation  such  as  a 
Laser  Doppler  Velocimeter  (LDV),  pulse  wire  anemometer,  or  split  film 
anemometer  may  have  been  better  suited  to  resolve  the  reverse  flow 
ambiguity  in  the  recirculation  region  behind  a  swept  backward-facing  step,  it 
was  felt  that  a  hot  wire  or  film  anemometer  was  capable  of  providing 
meaningful  data  within  the  scope  of  this  study.  Multi-component  probes 
were  not  used  for  several  reasons:  first,  none  were  available  and  they  were 
not  as  easy  to  use  as  a  single  component  probe;  second,  interference  between 
multiple  probes  being  spaced  too  close  together  could  cause  measurement 
problems;  third,  it  was  desirable  for  the  measurement  to  be  spacially  as  small 
as  possible  (i.e.,  point]  2md  since  the  size  of  many  multiple  component  probes 
were  considerably  larger,  the  single  component  probe  best  met  this 
requirement;  and  fourth,  the  flow  was  assumed  to  be  basically  two 
drmensionsd  (v  component  being  small),  thus,  within  this  assumption,  a 
single  probe  could  easily  make  the  necessary  measurements.  Because  of 
these  reasons  and  the  added  advantages  of  ruggedness,  abundant  supply  and 
relative  ease  of  use,  the  1210-20  Hot  Film  single  component  probe  was 
chosen.  Also,  since  the  intent  of  this  study  was  to  be  a  preliminary  attempt 
at  determining  the  mean  flowfield  behind  a  backward-facing  swept  step,  the 
disadvantages  of  using  a  hot  film  probe,  rather  than  a  hot  wire,  were  deemed 
acceptable.  The  drawback  of  using  a  hot  film  anemometer  is  a  lower 
frequency  response  but  since  only  mean  velocities  were  measured  this  was 
not  a  major  concern. 

Several  assiimptions  were  made  pert^ing  to  the  use  of  a  hot  film 
anemometer.  According  to  TSI  (TSI  Bulletin  TB18),  the  effects  of  heat 
transfer  from  the  probe  to  the  substrate  and/or  probe  supports  was  assumed 
to  be  negligible.  As  described  by  Lomas  (1986),  changes  in  humidity  and 
near  wall  effects  could  be  neglected.  The  closest  measurement  taken  to  the 
wall  was  0.020  inches  and  wall  effects  normally  do  not  become  appreciable 
imtil  a  distance  of  0.010  inches  from  the  wall.  Thus,  these  effects  were 
neglected.  Ambient  temperatme  changes  were  corrected  using  a  simple 
ratio  of  temperature  and  output  voltages  as  is  suggested  by  Spring  (1987) 
amd  TSI  TBS,  TB16  and  TB18.  This  correction  factor  will  be  presented 
when  the  data  reduction  procedure  is  discussed  in  chapter  4. 

Description  of  Data  Acquisition  Equipment 

A  Hewlett-Packard  data  acquisition  system  was  used  to  acquire  the 
calibration  and  velocity  profile  data  for  this  study.  This  system  consists  of 
the  following  hardware:  HP  3497 A  Controller,  which  switches  the  various 
data  channels  into  the  HP3456A  Digital  Integrating  Voltmeter  (it  in  turn, 
converts  the  hot  film  anemometer  bridge  output  voltages  from  amalog  to 
digital);  and  a  HP9826  Technical  Computer  system,  which  controls  the  data 
acquisition  system.  The  HP9826  is  programmed  in  HP  BASIC  and  contains 
1.5  mbytes  of  RAM. 
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Other  equipment  used  during  the  wind  tunnel  tests  included  a  0.25 
inch  diameter  pitot  static  tube,  which  was  used  to  measure  the  freestream 
velocity  for  calibrating  the  hot  film  anemometer.  The  pitot  static  tube  was 
connected  to  a  set  of  three  Magnehelic  Pressure  Gauges,  with  a  range  of  0.0 
to  0.25,  0.0  to  1.0  and  0.0  to  10.0  inches  of  water,  respectively.  The  ambient 
test  conditions  were  determined,  using  a  digital  thermocouple  to  measure  the 
test  section  air  temperature  and  a  barometer  to  measiire  atmospheric 
pressure.  The  ambient  temperature  and  pressure  were  required  to  calculate 
the  air  density  in  the  test  section. 
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Chapter  3 

EXPERIMENTAL  PROCEDURE 
Overview 


The  measTirement  of  the  flowfield  in  the  vicinity  of  a  backward-facing 
step  was  conducted  in  two  parts.  The  first  part  consisted  of  experimental 
flow  observations  and  velocity  measurements,  and  the  second  pzirt  was  the 
data  reduction  of  these  observations  and  measurements.  The  testing  and  the 
data  acquisition  procedure  are  presented  in  this  chapter.  Data  reduction  will 
be  discussed  in  Chapter  4. 

The  experimental  portion  of  this  study  was  also  aiccomplished  in  two 
parts,  the  first  part  being  a  qualitative  analysis,  using  oil  flow  visualization, 
of  the  flowfield  behind  the  three  backward-facing  step  configurations.  The 
second  part  of  the  test  used  a  single  component  hot  Him  anemometer  to 
provide  quantitative  measurements  of  the  velocity  profiles  at  different 
longitudinal  centerline  distances  both  upstream  and  downstream  of  the  step. 

Two  Reynolds  niimbers  definitions  usually  associated  with  the 
measurements  of  the  flowfield  over  a  backward-facing  step.  The  first 
Reynolds  number,  is  referenced  to  the  step  height,  Re]^,  and  the  other  to  the 
separation  location,  Re^ep.  For  this  study,  a  0.5  inch  step  height  yields  a 
value  for  Reii=2xl0^  and  with  a  distance  from  the  leading  edge  to  the  step 
of  27.0  inches  the  Reynolds  number  based  on  Lsq,  is  Resep=lxlO®.  These  two 
quantities  are  useful  when  comparing  results  from  different  studies. 

Oil  Flow  Visualization  Procedure 


The  oil  flow  experiments  were  conducted  to  determine  if  the  flow 
downstream  of  the  backward-facing  step  was  approximately  the  same  as 
presented  by  Selby  (1982).  The  oil  mixture  used  for  these  tests  was 
produced  by  combining  80.0  ml  of  100  centistoke  silicone  fluid  with  2.0  ml  of 
titanium  dioxide.  The  titanium  dioxide  does  not  dissolve  in  the  silicone  fluid 
but  remains  suspended.  The  mixture  was  poured  onto  the  plate  several 
inches  downstream  of  the  step.  The  tuimel  was  turned  on  and  the 
movement  of  the  fluid  observed.  The  results  from  the  oil  flows  downstream 
of  the  step  are  presented  in  Chapter  5. 

Oil  flow  observations  were  also  made  along  each  of  the  three  leading 
edges  of  the  three  swept  step  configurations.  Initially,  the  flow  over  the 
leading  edge  of  each  step  was  separated  and  the  downstream  velocity  profiles 
were  thus  2ifi'ected.  In  order  to  improve  the  flow,  each  step  configuration 
lezuling  edge  was  blxmted  to  a  0.25  inch  radius  and  hcuid  sanded  to  smooth 
the  transition  from  the  lower  plate  wedge,  around  the  radius,  and  then  onto 
the  upi>er  surface.  This  improved  the  flow  slightly  but  not  to  an  acceptable 
level.  Finally,  the  model  was  placed  at  a  slight  negative  angle-of-attack, 
approximately  2.0  degrees  nose  down.  This  greatly  improved  the  flow  over 
the  leading  edge  of  the  plate.  By  placing  the  plate  at  an  angle-of-attack  a 
small  favorable  pressure  gradient  was  generated.  After  excimining  the  free 
stream  velocities  at  severzd  positions  along  the  tunnel  centerline  little  effect 
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due  to  the  pressure  gradient  on  the  flow  could  be  detected.  A  grit  strip,  of 
#60  grit,  1.0  inch  wide,  was  placed  2.0  inches  downstream  of  each  leading 
edge  in  order  to  assme  transition  to  turbulent  flow  as  far  upstream  of  the 
step  as  possible. 


Hot  Film  Calibration  Procedure 


Calibration  of  the  probe  at  a  specific  ambient  temperature  took 
approximately  twenty  minutes  to  complete.  The  procedure  required  some 
manual  operations  such  as  probe  alignment  and  positioning,  and  tunnel 
velocity  ^justments,  but  the  data  were  acquired  and  stored  automatically. 
To  calibrate  the  hot  film  anemometer,  the  probe  was  positioned  half  way 
between  the  model  and  the  tunnel  ceiling,  approximately  10.0  inches  above 
the  wind  tunnel  model  surface.  The  tuimel  was  then  turned  on  and  set  at  70 
mph  (or  103  fps).  The  same  traversing  mechanism  used  for  the  velocity 
profile  tests  was  also  used  to  calibrate  the  probe. 

Hot  Film  Probe  Alignment 

Alignment  of  the  probe  was  accomplished  first  by  sight  positioning 
the  film  sensor  as  close  to  parallel  to  the  freestream  as  possible.  Then  the 
minimum  output  voltage  was  found  by  monitoring  the  bridge  output  voltage 
on  the  digital  voltmeter  while  making  fine  rotational  adjustments  about  the 
longitudinal  axis  of  the  probe.  The  minimum  output  voltage  corresponded 
to  the  film  sensor  being  aligned  parallel  to  the  flow.  When  the  minimum 
was  foimd,  the  yaw  angle  indicator  on  the'  traversing  adapter  was  set  to  90 
degrees.  As  a  check,  the  probe  was  then  rotated  plus  and  minus  10  degrees 
about  90  degrees  and  the  bridge  output  compared.  If  the  probe  was  properly 
aligned  the  two  readings  would  be  identical,  if  not  the  probe  alignment 
process  was  repeated.  With  the  alignment  completed,  the  probe  was  then 
rotated  so  that  the  yaw  angle  indicator  read  0  degrees  and  the  film  sensor 
was  perpendicular  to  the  flow.  The  next  step  was  to  record  the  effect  of 
velocity  on  the  bridge  output  voltage  of  the  hot  film  probe. 

Velocity  Calibration 


The  velocity  calibration  begcui  with  the  input  of  initial  conditions  into 
the  data  acquisition  program.  These  inputs  included:  test  and  run  numbers, 
date  of  the  test,  comment  about  the  test,  probe  serial  number,  test  section 
temperature,  ambient  pressure,  calibration  velocity  range,  bridge  reference 
voltage,  auid  ntimber  of  readings  per  sample.  The  calibration  velocity  range 
was  70  mph  (103  fps)  to  7  mph  (10  fps).  The  wind  tunnel  manometer  board 
was  used  only  as  a  reference  to  set  the  approximate  test  section  velocity. 

The  actual  velocity  was  determined  by  manually  reauiing  the  appropriate 
Magnehelic  water  gauge  which  was  connected  to  the  pitot  tube  installed  in 
the  test  section.  Starting  at  70  mph,  the  data  acquisition  system  read  the 
bridge  output  voltage  and  stored  the  average  value  of  50  readings.  The 
voltmeter  was  set  for  a  reading  rate  of  50  readings  per  second;  thus,  it  took 
approximately  1.0  second  to  complete  one  sample.  The  test  section  velocity 
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was  decreased  by  5  mph  and  the  process  repeated  until,  a  velocity  of 
approximately  5  mph  was  reached. 

Yaw  Calibration 


The  yaw  calibration  of  the  probe  was  the  next  step  in  the  calibration 
process.  The  probe  was  yawed  by  rotating  the  probe  about  its  longitudinal 
axis.  In  order  for  a  single  component  probe  to  resolve  the  direction  of  the 
total  velocity  vector,  the  effects  of  yaw  angle  on  the  probe  output  must  be 
known.  The  reason  for  this  yaw  effect  is  that  as  the  yaw  angle  increases,  the 
component  of  the  velocity  perpendicular  to  the  film  sensor  decreases.  Thus, 
the  effective  cooling  of  the  sensor  and  the  bridge  output  voltage  decreases. 
Generedly  the  yaw  effect  on  a  hot  film  or  wire  follows  a  cosine  law,  defined  as 
follows: 


Ueff  =  U  cos(q!)  (3.1) 

but  according  to  Lomas  (1986)  thr  is  only  an  approximation  if  the  probe’s 
aspect  ratio  is  less  than  600.  The  probe  aspect  ratio  is  defined  as  the  atio  of 
the  length  to  diameter  o^  the  film  sensor.  For  this  study  the  hot  film  probe 
had  an  aspect  ratio  of  about  20.  Since  the  aspect  ratio  was  so  low,  probe  sup¬ 
port  interference  could  also  affect  the  accuracy  of  the  cosine  law.  Thus,  the 
yaw  calibration  of  the  hot  film  probe  was  necessary  for  these  tests. 

The  yaw  calibration  was  conducted  at  a  velocity  of  50  mph  (73,3  fps), 
which  was  aJso  the  freestream  velocity  used  for  all  the  velocity  profile  surveys. 
It  should  be  noted  that  the  effect  of  yaw  on  the  non-dimension^  probe  output 
is  Independent  of  the  freestream  velocity,  as  stated  by  Olivari  (1978),  and  was 
ejperimentally  verified  during  this  study.  The  yaw  c^ibration  was  performed 
with  the  probe  positioned  at  the  same  height  above  the  step  as  was  used  for 
the  velocity  calibration.  The  probe  was  rotated  to  a  -60  degree  yaw  angle  and 
the  average  of  50  readings  of  the  bridge  output  was  recorded  and  stored.  The 
yaw  angle  was  incremented  10  degrees  and  the  process  repeated  until  a  +60 
degree  yaw  zuigle  was  achieved. 

All  calibration  data  were  stored  on  a  flexible  disc  for  later  reduction 
and  analysis.  Whenever  the  test  section  temperature  for  a  velocity  profile  test 
was  more  than  three  degrees  Fahrenheit  from  euiy  of  the  calibration  tempera¬ 
tures,  a  new  calibration  was  performed.  The  calibration  results  and  data 
reduction  technique  are  discussed  in  the  Chapter  4. 

Velocity  Profile  Test  Procedure 

The  test  procedured  for  performing  the  velocity  profile  sxirvey  was 
similzir  to  the  calibration  test.  First,  the  initial  conditions  for  the  data  acquisi¬ 
tion  program  were  entered.  These  included:  a  test  and  run  number,  test  sec¬ 
tion  temperature  and  cunbient  pressure,  sweep  angle  of  the  step,  position  along 
the  model  centerline  with  respect  to  the  step  where  the  profile  would  be  sur¬ 
veyed,  and  the  probe  rotation  or  yaw  angle,  which  was  arbitrary  chosen  to  be 
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15  degrees.  The  traversing  mechainism  was  positioned  and  locked  into  place 
above  the  point  where  a  velocity  profile  was  to  be  measured. 

The  next  step  was  to  establish  a  reference  height  between  the  surface 
of  the  plate  and  the  tip  of  the  hot  film  probe.  It  was  discovered  that  because 
of  the  negative  pressure  in  the  test  section  during  tunnel  operations,  the  tun¬ 
nel  ceiling  and  floor  would  move  relative  to  each  other.  This  made  it  difficult 
to  establish  a  reference  height  between  the  probe  tip  and  the  surface  of  the 
model.  It  was  also  discovered  that  other  unknown  factors,  possibly  tempera- 
tiire,  affected  the  amoimt  of  displacement.  After  extensive  testing  the  follow¬ 
ing  procedure  was  developed  to  provide  the  necessary  accuracy  in  the  height 
measurement.  This  technique  established  the  height  of  the  probe  above  the 
surface  of  the  plate  to  plus  or  minus  0.002  inches.  The  technique  consisted  of 
replacing  the  probe  in  the  probe  holder  with  a  measurement  standard,  a  modi¬ 
fied  hot  film  probe.  The  exact  difference  in  length  between  the  measurement 
standard  and  the  hot  film  probe  was  determined  using  cin  optical  comparator. 
The  wind  tunnel  was  turned  on  and  the  velocity  adjusted  to  the  test  speed  of 
50  mph.  After  allowing  the  tunnel  to  run  for  several  minutes  the  standard 
was  lowered  to  the  surface  of  the  model  «md  the  height  indicator  on  the 
traversing  mechanism  was  zeroed.  The  standzird  was  raised  and  lowered  imtil 
the  height  indicator  read  zero  repeatedly.  The  standard  was  then  offset  from 
the  plate  a  distance  equal  to  the  difference  in  length  of  the  standard  to  the 
hot  film  probe.  The  height  indicator  was  then  reset  to  zero.  This  position 
then  indicated  the  point  at  which  the  probe  tip  would  touch  the  plate.  With 
the  height  reference  set,  the  tunnel  was  turned  off  and  the  measurement  stan¬ 
dard  replaced  with  the  hot  film  probe. 

Next  the  probe  was  positioned  10.000  inches  above  the  surface  of  the 
plate.  At  this  position  the  probe  was  aligned  with  the  freestream  velocity 
direction  using  the  same  alignment  procedure  as  discussed  above  for  the  velo¬ 
city  calibration  of  the  probe.  It  was  assumed  that  at  10.000  inches  above  the 
plate  the  freestream  velocity  was  parallel  to  the  sides  of  the  tunnel.  The  velo¬ 
city  surveys  consisted  of  three  measurements  taken  at  ezurh  vertical  position 
above  the  plate.  Only  two  measurements  were  necessary  to  resolve  the  mag¬ 
nitude  and  direction  of  the  total  velocity,  and  the  third  measTirement 
increased  the  accuracy  of  this  procedme.  The  first  measurement  was  taken 
with  the  yaw  angle  indicator  set  at  zero  degrees.  The  average  of  50  readings 
was  recorded  and  stored.  The  probe  was  then  rotated  plus  and  minus  15 
degrees  about  0  degrees,  with  the  average  of  the  50  reading  being  read  and 
stored  for  each  rotation.  The  same  process  was  repeated  for  each  vertical  sur¬ 
vey  position  above  the  plate.  Vertical  positioning  varied  depending  on 
whether  the  profile  was  being  taken  upstream  of  the  step,  in  the  recirculation 
region,  or  downstream  of  the  recirculation  region.  Most  profiles  began  at 
10.000  inches  and  ended  at  0.020  inches  above  the  step.  The  lower  limit  of 
0.020  inches  was  chosen  in  order  to  avoid  breaking  probes  by  inadvertently 
driving  them  into  the  plate.  The  same  procedure  was  then  repeated  for  vari¬ 
ous  longitudinal  positions  upstream  and  downstream  for  e<ich  of  the  three  step 
configurations. 
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Flow  Reversal  Ambiguity 


For  this  first  attempt  at  measuring  the  flowfield  in  the  vicinity  of  a 
swept  backward-facing  step,  no  provisions  were  made  to  experimentcilly 
resolve  the  flow  reversal  ambiguity  within  the  recirculation  region.  Instead, 
the  point  of  flow  reversal  was  determined  during  the  data  reduction  pro¬ 
cedure,  the  assumption  being,  that  as  the  stagnation  line  in  the  recirculation 
region  was  approached  from  the  freestream  side,  the  total  velocity  would 
decrease  to  a  minimiim.  Then,  on  the  model  side  of  the  stagnation  line,  the 
total  velocity  would  increase  slightly  before  going  to  zero  again  at  the  surface 
of  the  model.  This  approach  was  only  partially  successful.  Two  simple  exper¬ 
imental  methods  may  prove  useful  in  futme  testing  for  resolving  the  flow 
reverssil  2Lmbiguity.  One  method  would  be  to  use  a  small  rotatable  pitot  probe 
to  find  the  direction  of  maximum  pressure  and  thus  the  flow  direction.  The 
other  method  described  by  Naverbury  (1969)  uses  a  hot  wire/film  with  a 
shield  around  it.  The  shield  causes  the  probe  to  be  sensitive  to  the  flow  in 
only  one  direction.  Then,  similar  to  the  pitot  probe,  by  rotating  the  hot 
wire/film  probe  to  find  the  maximum  bridge  output  voltage,  the  flow  direction 
can  be  determined.  However,  for  these  preliminary  tests  these  techniques 
were  not  attempted. 


Chapter  4 

DATA  REDUCTION  PROCEDURE 
Overview 


This  chapter  presents  a  discussion  of  the  procedure  used  in  this  study 
to  reduce  the  velocity  profile  measurements.  The  calibration  procedure, 
including  the  velocity  and  yaw  efiects  on  the  hot  film  probe  are  discussed 
first.  Next  the  data  reduction  procedure  to  determine  the  total  velocity  and 
flow  angle  from  the  three  measurements  at  each  survey  position  is  explained. 
The  data  analysis  procedure  for  determining  the  flow  reversal  point,  the  u 
and  w  velocity  components  and  the  u  component  boundary-layer  thickness  is 
reviewed.  Typical  results  from  each  procedure  are  presented.  Finally,  the 
format  for  the  tabulated  and  graphical  presentation  of  the  results  from  the 
above  procedures  are  explained.  The  intermediate  results  from  the  data 
reduction  procedure  have  not  been  included  for  the  sake  of  brevity.  The 
intermediate  results  can  be  obtained  upon  request. 

Calibration  Data  Reduction  Procedme 


Velocity  Calibration  Reduction  Procedure 

The  procedure  followed  to  calibrate  the  hot  film  probe  for  the  effect 
of  velocity  on  bridge  output  was  relatively  straightforwiird.  The  acquisition 
of  the  calibration  data  was  described  in  Chapter  3.  With  the  film  sensor 
aligned  perpendicular  to  the  freestream  the  bridge  output  from  the 
anemometer  was  recorded  along  with  the  true  test  section  velocity,  which 
was  determined  from  a  pitot  static  probe.  The  velocity  calibration  was 
calculated  using  two  different  methods. 

Fourth  Order  Polynomial  Curve  Fit.  The  first  approach  fitted  the 
raw  velocity  and  output  voltage  data  to  a  foxirth  degree  polynomial,  using  a 
least  square  routine.  This  approach  is  similar  to  the  way  the  1052  Linearizer 
component  of  the  TSI  Hot  Wire  Anemometer  System  functions.  The 
difference  is  that  the  Linearizer  fits  a  fourth  degree  polynomial  to  the  bridge 
output  electronically  where  this  procedure  fits  the  polynomial  nximerically. 

A  typical  plot  of  the  raw  calibration  data  and  resulting  curve  fit  is  shown  in 
Figure  4.1. 

Linear  Ci^e  Fit.  The  second  method  of  determining  a  relationship 
between  the  bridge  output  and  the  velocity  was  by  using  King’s  law  zmd  a 
least  squares  linear  fit  of  the  data.  King’s  law  describes  the  heat  transfer 
from  a  cylinder  of  infinite  length.  The  law  for  a  hot  film  anemometer  is 
often  written  as 


E^  =  A  +  BU“ 


(4.1) 


where  E  is  the  hot  wire  bridge  output  voltage,  U  is  the  velocity,  A  and  B  are 
constants  and  n  varies  between  0.45  and  0.5.  For  this  study,  n  =  0.5  provided 
the  best  results.  A  and  B  were  determined  using  a  linecir  least  squcires  fit  of 
the  calibration  data.  Figure  4.2  compares  the  experimental  data  to  linear 
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(VELOCITY)-^.  5  C-Fps^.5)  BRIDGE  OUTPUT 


Figure  4.1  Raw  Hot  Film  Calibration  Data  and  Fourth 
Degree  Polynomial  Curve  Fit 


Figure  4.2  King's  Law  Applied  to  Hot  Film  Calibration 
Data  and  Compared  to  Linear  Curve  Fit 
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curve  fitted  data  for  a  typical  calibration  test.  It  should  be  noted  that  the 
measurement  taken  at  zero  velocity  was  not  used  in  the  least  squares  calcula¬ 
tion  because  of  the  potential  for  buoyancy  effects.  Better  agreement  was 
achieved  between  the  true  and  calculated  velocities  over  the  entire  calibration 
range  by  neglecting  the  zero  velocity  measurement. 

Comparison  of  Polynomi^  and  Linear  Curve  Fits.  Typical  calibration 
data  is  presented  in  Table  4.1.  The  important  information  in  this  table  is  the 
percent  difference  between  the  velocity  measured  with  the  pitot  static  probe, 
and  the  velocities  calculated  by  the  two  curve  fitted  equations:  linear  and 
fourth  order  polynomial.  Overall  the  agreement  was  very  good,  with  percent 
differences  generally  around  2%  or  less.  At  the  lower  velocities,  the  percent 
difference  was  slightly  higher  with  the  highest  being  5%,  but  this  was  to  be 
expected.  The  increase  in  error  at  the  low  velocities,  below  10  mph,  was  due 
to  the  reduced  sensitivity  of  the  water  gauges  at  such  low  pressures.  The 
uncertainty  of  the  measxirements  is  disciissed  in  the  last  section  of  this 
chapter.  All  velocities  were  calculated  with  the  polynomial  fit,  since  it  was 
slightly  more  accurate  than  the  linear  curve  fit. 

Yaw  Calibration 


Eff^t  of  Probe  Yaw  on  Velocity  Measurem^t.  Typical  yaw  charac- 
teristics  of  a  hot  film  probe  are  listed  in  Table  4.2.  The  first  step  in  determin¬ 
ing  the  effects  of  yaw  angle  on  the  probe’s  output  was  to  calculate  the  veloci¬ 
ties  corresponding  to  the  yawed  voltage  readings,  obtained  during  calibration, 
using  the  fourth  order  polynomial  from  the  velocity  calibration.  Next,  the  yaw 
data  were  non-dimensionaJized  by  dividing  the  velocities  by  the  maximum 
velocity  which  occurred  at  0  degrees.  Then,  the  yaw  angle  and  the  non- 
dimensional  velocity  were  curve  fitted  using  a  second  order,  least  squares  tech¬ 
nique.  Fig\ire  4.3  compares  the  experimental  data  with  the  resulting  yaw  cali¬ 
bration  curve.  As  can  be  seen  from  the  plot  and  the  tabulated  percent  differ¬ 
ences  the  agreement  is  quite  good,  2%  error  or  less.  From  this  c\irve  the  total 
velocity  could  be  determined  once  the  flow  direction  was  known. 

Velocity  Profile  Data  Reduction 

A  procedure  for  determining  the  total  flow  velocity  and  direction  in  an 
iinknown,  two-dimensional  flowfield  was  described  by  Olivari  (1978).  In  this 
study,  the  two  dimensions  of  interest  were  the  components  perpendicular  (u) 
and  parallel  (w)  to  the  step  face.  The  v  component  of  the  velocity  vector  was 
assumed  to  be  negligible.  An  overview  of  the  procedure  will  be  given  first,  fol¬ 
lowed  by  an  explanation  of  how  it  was  used  to  resolve  the  total  velocity  and 
direction  from  the  velocity  profile  measurements. 
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Table  4.1  Typical  Hot  Wire  Calibration  Results 


HOT  UIRE  CALIBRATION  RESULTS 
DATE  :  1Z/Z0/88 
TEST  NO. -RUN  NO.:  IZI0-Z 
COMMENT:  CALIBRATION  RUN 
PROBE  SERIAL  NO.:  a5Z34 
CALIBRATION  TEMPATURE:  58  den  F 
PRESSURE:  764.4  nn-hg 
SAMPLES/ READING:  50 
TEST  FREE  STREAM  VELOCITY:  50  fps 
CALIBRATION  VELOCITY  RANGE:  5  TO  70  noh 
VOLTAGE:  5.004ZS  volts 
YAW  CALIBRATION  ANGULAR  STEP  10  deg 


N 

IND  VEL 

PRESSURE 

nph 

inch-HZO 

1 

0.0 

0.000 

Z 

70.0 

Z.600 

3 

65.0 

Z.Z00 

4 

60.0 

1.300 

5 

55.0 

1.600 

6 

50.0 

1.350 

7 

45.0 

1  .  100 

8 

40.0  • 

.870 

3 

35.0 

.695 

10 

30.0 

.5Z0 

1  1 

ZS.0 

.330 

12 

Z0.0 

.Z65 

13 

15.0 

.  148 

14 

10.0 

.091 

15 

8.0 

.053 

IS 

5.0 

.0ZZ 

ACTUAL 

VELOCITY 

BRIDGE 

BRIDGE 
VEL  LIN 

fps 

volts 

fps 

0.00 

3. 1ZZ 

.81 

106.  13. 

6.073 

107.86 

37.68 

5.974 

33.33 

30.77 

5.878 

31.61 

83.30 

5.770 

83.38 

7S.51 

5.66Z 

75.68 

69.07 

5.546 

67.39 

61 .4Z 

5.413 

53.86 

54.90 

5.303 

53.  BZ 

47.49 

5. 17Z 

46. 86 

41 . 13 

5.037 

40.50 

3J.90 

4.865 

33.27 

ZS.33 

4.673 

ZB. 49 

13.87 

4.430 

Z0.61 

16.00 

4.344 

16.73 

9.77 

3., 936 

9.46 

DELTA 

BRIDGE 

DELTA 

2 

VEL  POLY 
fps 

2 

0.0 

.05 

0.0 

t:b 

T05’.4T — 

-  .7 

1  .7 

38.  12 

.5 

.3 

91.29 

.6 

.  1 

83.76 

.  b 

- 1 . 1 

76.50 

-0.0 

-1.6 

63.03 

-0.0 

-Z.5 

60.33 

-.8 

-2.3 

54.58 

-.5 

-1.3 

47.58 

.  2 

-1.5 

40.93 

-.5 

-1.9 

33.31 

-1.7 

4.5 

26.  13 

3.4 

3.7 

20.  1  3 

1  .  3 

4.6 

16.22 

1.4 

-3.  1 

9.28 

-5.0 

Table  4.2  Typical  Yaw  Calibration  Results 


YAW  CHARACTERISTICS,  based  on  bridge  output 

ANGLE 

YAW 

YAU  VEL 

UE/U 

UE/U 

X 

deg 

60 

volts 

5.077 

fps 

4Z.8Z 

actual 

.558 

poly 

.556 

-.5 

50 

5.ZB4 

52. 46 

.684 

.630 

.3 

40 

5.4Z1 

61 .40 

.300 

.800 

-.  1 

30 

5.544 

68.87 

.898 

.885 

-1.4 

20 

5.585 

71.48 

.932 

.346 

1 . 6 

10 

5.648 

75.56 

.385 

.383 

-.2 

0 

5.665 

76. 7Z 

1.000 

.996 

-.4 

-10 

5.651 

75.79 

.388 

.384 

-.4 

-20 

5.585 

71  .52 

.93Z 

.348 

I  .7 

-30 

5.535 

68.35 

.891 

.887 

-.4 

-40 

5.447 

SZ.31 

.820 

.803 

-2.  1 

-50 

5.264 

52.47 

.684 

.634 

1 .4 

-60 

5.080 

42.35 

.560 

.560 

.  1 

NOTE: 

YAU_VEL 

's  are  calculated 

using  a 

polynoninal 
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Second  Order  Polynomial  Curve  Fit 


Overview  of  Data  Reduction  Procedure 


By  using  the  yaw  calibration  curve,  which  represents  the  angular  sensi¬ 
tivity  of  the  hot  film  probe,  the  magnitude  and  direction  of  the  velocity  in  an 
unknown  two-dimensionzd  flowfield  can  be  found.  At  ecich  test  point  in  the 
velocity  profiles,  three  measurements  are  made,  a,  a  +  9,  and  a  —  9,  where 
a  is  the  unknown  flow  direction  relative  to  some  reference,  suid  d  is  eui  arbi¬ 
trarily  chosen  angle.  The  only  requirement  for  this  procedure  is  that  the 
angles  a,  a  +  9  and  a  —  0  be  less  than  60  degrees,  because  the  calibration 
curve  is  only  valid  between  -60  and  +60  degrees.  In  this  study,  the  reference 
was  chosen  to  be  the  free  stream  flow  direction.  The  following  three  equations 
represent  the  three  measxirements: 


U 

Ueff, 

u 

Ueff. 

U 


(4.2) 

II 

+ 

(4.3) 

=  l{a-e) 

(4.4) 

The  function  f  is  the  yaw  calibration  curve  determined  from  the  reduction  of 
the  calibration  data.  The  next  step  is  to  ratio  equation  (4.2)  to  (4.3)  and  (4.4) 
to  (4.2).  These  ratios  yield  the  following: 


Urfr./U  f{a  +  «)  ^  ' 


(4.5) 


f(a) 


G{a) 


(4.6) 


Equations  (4.5)  and  (4.6)  allowed  the  flow  direction,  a ,  to  be  determined 
because  the  three  hot  film  measurements  and  9  are  known.  Actually,  two 
separate  a ’s  are  found  2uid  averaged  to  give  the  final  flow  direction  at  that 
particxilar  point.  Knowing  the  average  a ,  it  is  then  substituted  into  (4.2)  to 
determine  the  magnitude  of  the  total  velocity.  If  the  requirement  for  a , 
a  9 ,  a  —  9  being  less  than  60  degrees,  is  not  met,  then  the  calculated  a  is 
only  an  estimate,  which  can  be  used  to  reposition  the  probe  to  find  a  more 
ar.curate  result. 
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Example  of  Data  Reduction  Procedure 


The  first  step  in  reducing  the  2uiemometer  voltage  measurements,  to 
velocity  and  flow  direction  is  to  compensate  for  temperature  differences 
between  the  cedibration  and  velocity  profile  tests.  This  compensation  is  neces¬ 
sary  because  of  the  way  in  which  a  constant  temperature  hot  film  anemometer 
operates.  The  film  sensor  makes  up  one  leg  of  a  Wheatstone  Bridge  and  is 
typically  heated  to  almost  500  degrees  Fahrenheit.  As  a  fluid,  such  as  air, 
flows  over  the  sensor,  the  sensor  is  cooled  and  the  bridge  becomes  unbalanced. 
The  anemometer  balances  the  bridge  to  maintain  a  constant  sensor  tempera¬ 
ture  by  adjusting  the  voltage  across  the  sensor,  thus  the  name,  constant  tem¬ 
perature  anemometer.  If  the  sensor  is  calibrated  at  one  temperature  and  then 
tested  at  another,  the  anemometer  will  adjust  for  the  temperature  difference 
by  rzusing  or  lowering  the  bridge  voltage.  Since  the  bridge  voltage  and  fluid 
velocity  are  related  by  the  calibration  equation  which  was  determined  at  a 
specific  calibration  temperature,  any  temperature  difference  will  cause  the 
anemometer  to  indicate  an  incorrect  fluid  velocity. 


There  are  several  ways  to  correct  for  temperature  variations  between 
calibrated  and  test  temperatures  when  using  a  hot  film  anemometer.  Most 
are  complicated  and  require  several  assumptions,  which  leave  the  results  ques¬ 
tionable.  A  simple  method  suggested  by  Spring  (1987)  and  TSI  (TB18), 
adjusts  for  the  temperature  difference  by  a  single  correction  factor. ,  The 
correction  factor  is  defined  as: 


Ecorr^  _ 

e2  ■  t,  -  t 


(4.7) 


where  Econ  is  the  co.rrected  bridge  output  voltage,  E  is  the  voltage  from  the 
velocity  profile  measurements,  Tg  is  the  sensor  temperature  supplied  by  the 
manufacture  (482  degrees  Fahrenheit),  T^ai  is  the  fluid  temperature  during 
calibration,  and  T  is  the  temperature  of  fluid  during  the  velocity  profile  tests. 
This  correction  factor  was  used  for  temperature  differences  up  to  3  degrees 
Fahrenheit.  If  the  difference  was  greater,  a  new  calibration  was  performed  at 
the  new  temperature. 

In  order  to  reduce  the  velocity  profile  data,  an  approximation  to  the 
functions  F( a),  equation  (4.5),  and  G(  a),  equation  (4.6),  had  to  be  deter¬ 
mined.  F(  a)  was  found  by  generating  a  locus  of  data  using  the  yaw  calibra¬ 
tion  curve.  Starting  with  an  a  equal  to  -60  degrees  and  equation  (4.2),  the 
velocity  ratio  corresponding  to  that  particular  adpha  was  calculated.  Next,  the 
corresponding  velocity  ratio  for  equation  U.3)  was  found  using  a  0  equal  to  15 
degrees.  Then  employing  equation  (4.5),  the  value  of  F  was  determined  for  an 
a  of  -60  degrees.  The  process  was  repeated  incrementing  a  by  -f-5  degrees 
until  an  a  of  4-45  degrees  was  achieved.  The  upper  limit  of  4-45  degrees  was 
required  in  order  that  9  +  a  did  not  exceed  4-60  degrees  because  the  yaw  cali¬ 
bration  curve  was  only  valid  between  -60  and  4-60  degrees.  With  the  locus  of 
points  calculated,  F(  a)  was  approximated  by  a  fourth  order,  least  squares 
polynomial.  A  slight  modification  to  the  above  procedure  was  necessary  in 
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order  to  approximate  G(  a).  The  initial  and  final  a  was  -45  degrees  and  +60 
degrees,  respectively.  B  was  still  15  degrees,  but  instead  of  equation  (4.5), 
equation  (4.6)  was  used.  Figure  4.4  compares  the  locvis  of  points  and  the 
fourth  order  polynomial  curve  fits  for  the  function  F(  a)  and  G(  a).  The 
agreement  between  the  data  and  curve  fits  were  very  good,  with  a  maximum 
error  of  less  than  1%. 

With  the  approximation  ofF(o£),G(a)  and  the  three  survey  measure¬ 
ments  at  a ,  a  +  ^  and  a  —  0,  the  flow  direction  ( a)  was  calculated.  Know¬ 
ing  a  and  the  effective  velocity  measured  at  a ,  and  using  the  yaw  czdibration 
equation  (4.2),  the  magnitude  of  the  total  velocity  was  determined. 

Data  Reduction  Procedure 


The  data  reduction  procedure  described  in  this  section  did  not  provide 
satisfactory  results  in  the  region  of  x/h=0.0  to  1.0  and  from  y/h=0  to  1.0. 

This  was  due  to  the  total  velocity  being  very  low  and  out  of  range  of  the  cali¬ 
bration  data.  Also,  the  flow  reversal  could  not  be  resolved  in  this  region. 

Flow  Reversal  Determination.  Once  the  magnitude  of  the  total  velocity 
and  the  flow  direction  were  calculated,  other  quantities  necessary  for  the 
analysis  of  the  flowfield  in  the  vicinity  of  the  backward-facing  swept  step  were 
determined.  These  quantities  included  the  u  and  w  velocity  components  and 
the  boundary-layer  thickness.  Although  the  velocity  profiles  were  measured 
to  a  height  of  10.000  inches  (20  step  heights)  above  the  surface  of  the  model, 
only  data  between  2.000  inches  (4  step  heignts)  and  the  surface  were  used  in 
the  analysis.  The  reason  for  this  is  that  the  upstream  and  downstream 
boundary-layers  are  much  smaller  than  2.0  inches  and  the  boundary-layer  and 
separated  flow  region  are  the  primary  focus  of  this  work. 

After  modifying  the  data  acquisition  program  to  assist  with  the  data 
analysis  part  of  the  study,  the  first  step  was  to  place  the  upper  limit  of  the 
profile  at  2.000  inches.  Next,  the  total  velocity  profile  was  displayed  on  the 
computer  CRT  to  determine  if  a  flow  reversal  was  present.  As  explained  in 
chapter  3,  the  flow  reversal  was  identified  by  examining  the  total  velocity  pro¬ 
file  and  looking  for  a  minimxim  followed  by  a  slight  increase  in  the  total  mag¬ 
nitude  of  the  velocity.  Figure  4.5  shows  the  toti  velocity  for  the  0  degree 
swept  step,  x=l  inch  downstream  of  the  step.  The  flow  reversal  occtirs 
approximately  0.25  inches  above  the  surface.  If  a  particular  profile  was  not  in 
the  recirculation  region  or  a  clear  flow  reversal  was  not  evident,  the  total  velo¬ 
city  profile  was  accepted  as  displayed  and  the  data  analysis  continued.  If  a 
flow  reversal  was  present,  as  is  shown  in  Figure  4.5,  all  data  points  below  the 
stagnation  line  were  modified.  This  modification  included  reversing  the  sign  of 
both  the  total  velocity  and  the  flow  direction.  Generally,  at  the  point  of  flow 
reversal,  the  total  velocities  were  close  to  or  below  the  lower  calibration  limit 
of  the  hot  film  probe,  so  that  the  accuracy  of  these  results  were  questionable. 
In  some  cases,  it  was  necessary  to  drop  one  or  two  data  points  aroxmd  the  flow 
reversal  point  because  of  the  lack  of  reliability  of  the  calibration  curve  at  low 
velocities.  Figure  4.6  is  the  profile  of  Figure  4.5  modified  for  flow  reversal  and 
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Figure  4,6  Total  Velocity  Profile  With  Flow 

Reversal.  A  =  0  Degrees  and  x/h  =  2.0 
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two  data  points,  corresponding  to  probe  heights  of  0.25  and  0.20  inches,  being 
deleted  from  the  plot.  It  should  be  noted  that  all  data  points  are  listed  in  the 
tabulated  data  with  the  points  that  axe  not  used  in  the  analysis  indicated  at 
the  bottom  of  the  table. 

Edge  u  Velocity  Component.  With  the  flow  reversal  determined,  a 
graph  of  tne  u  velocity  component  was  displayed.  This  was  necessary  for 
determining  the  edge  free  stream  u  velocity  component.  The  edge  velocity 
was  determined  by  averaging  the  velocities  of  all  data  points  between  y=2.0 
inches  and  below  to  the  first  data  point  that  appeared  to  be  in  the  botmdary- 
layer.  A  typical  plot  of  the  u  velocity  component  versus  probe  height  is 
shown  in  Figure  4.7  with  the  edge  u  component,  u«,  indicated  by  the  dashed 
line.  The  edge  u  velocity  component  was  needed  for  several  reasons:  first,  the 
u  and  w  components  were  both  non-dimensionalized  by  u^ ;  secondly,  the  cut¬ 
off  for  determining  u«  was  also  used  to  calculate  zm  edge  total  velocity,  Ug,  for 
the  purpose  of  non-dimension<ilizing  the  total  velocity  results;  and  finally,  Ug 
was  needed  to  approximate  the  boimdary-layer  thickness. 

Approximate  Boundary-Layer  Thickness.  The  final  quantity  deter¬ 
mined  was  an  approximate  boimd^-layer  thickness  for  the  u  component 
velocity  profile.  By  employing  a  power  law  and  using  a  In-ln  plot  an  estimate 
of  the  boundary-layer  thickness  was  made.  The  computer  displayed  a  plot  of 
ln{u)  verses  ln(y) ,  this  is  shown  in  Figure  4.8.  The  edge  velocity  is  indicated 
by  a  horizontal  line  through  the  data  points  that  were  used  to  determine  the 
average.  All  remaining  velocity  data,  below  the  last  data  point  used  to  deter¬ 
mine  u«  and  above  the  data  point  where  a  possible  flow  reversal  took  place, 
were  used  in  a  linear  least  squares  fit.  The  linear  fit  is  also  indicated  on  the  In 
plot.  Two  criteria  used  for  selecting  the  best  linear  fit  were:  first,  how  well 
did  the  curve  fit  the  available  data  and  second,  where  did  the  linear  fit  cross 
the  horizontal  line  representing  the  edge  velocity.  If  the  fit  was  unsatisfeictory 
or  it  did  not  intersect  the  constant  ln(u,)  line  properly,  the  number  of  data 
points  in  the  linear  fit  was  adjusted.  The  criterion  for  the  proper  intersection 
of  the  two  lines  was  that  the  linear  ciirve  fit  had  to  intersect  the  horizontal  u^ 
line  near  the  last  data  point  used  to  determine  u^.  Once  the  graph  of  the  two 
lines  was  acceptable,  the  intersection  of  the  two  lines  determined  the  value  of 
the  ln(  5).  Figure  4.9  shows  an  example  of  the  final  results  from  this  pro¬ 
cedure.  With  this  final  piece  of  information,  the  results  were  tabulated  and 
plotted. 

Format  for  Presentation  of  Results 


Tabulated  Results.  The  results  from  the  velocity  profile  measurements 
are  presented  in  both  tabulated  and  graphical  formats.  Typical  tabulated 
results  are  listed  in  Table  4.3  for  a  step  sweep  2mgle  of  15  degrees  at  a  distance 
of  x=l  inch  downstream  of  the  step.  At  the  top  of  each  table,  the  test  condi¬ 
tions  are  listed,  along  with  the  edge  or  average  total  velocity,  the  edge  or  aver¬ 
age  u  component  of  the  total  velocity  and  the  boimdary  layer  thickness  of  the 
u  component.  All  results  were  temperature  compensated  in  order  to  account 
for  any  temperature  difference  between  the  calibration  test  and  the  velocity 
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Table  4.3  Typical  Velocity  Profile  Results 


VELOCITY  PROFILE  RESULTS 

DfiTE  :  i 2/ 14/83 

TEST  NO. -RUN  NO.:  1203-3 

COKHENT: VELOCITY  PROFILE  fiLPHn-i  DEB  5B 

TUNNEL  VELOCITY  FROM  PITOT  PROBE:  7S.75 

STEP  ANGLE:  15  degrees 

DISTANCE  FROM  STEP:  1  inches 

CALIBRATION  TEST  AND  RUN  NO.:  CAL12032 

CALIBRATION  TEnPATURE:  44  deg  F 

TEST  TEMPERATURE:  45  deg  F 

TEST  PRESSURE:  7BB.5  nn-Hg 

TEST  DATA  TAKEN  AT  -r  AND  -  VALUES  OF  15 

EDGE  TOTAL  VELOCITY':  73.S32  fps 

EDGE  u  COMPONENT  OF  VELOCITY':  74.564  fps 

BOUNDARY  LAYER  THICKNESS:  .552  inches 


DATA  TEMPERATURE  COMPENSATED 
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profile  test.  The  tabulated  results  include  the  following:  y,  distance  from  the 
probe  tip  to  the  model  surface  in  inches;  U/Ug,  total  velocity  non- 
dimensionalized  by  the  edge  total  velocity;  a ,  flow  direction  with  respect  to 
the  freestream;  o:  +  A,  positive  sense  in  the  clockwise  direction;  u/ug,  the 
component  of  the  total  velocity  normal  to  the  step  face  non-dimensionalized 
by  the  u  component  edge  velocity,  positive  sense  is  in  freestream  direction; 
w/ug,  the  component  of  the  total  velocity  parallel  to  the  step  face  non- 
dimensionalized  by  the  u  component  edge  velocity,  positive  sense  is  right  to 
left  when  looking  downstream.  The  total  velocity  was  determined  from  the 
fourth  order  polynomial  which  was  more  accurate  than  the  linear  fit  of  King’s 
law. 


Plotted  Results.  The  data  reduction  and  analysis  program  plotted  the 
results  on  seven  different  graphs.  An  example  of  these  plots  for  the  15  degree 
swept  step,  at  x=l  inch  downstream  of  the  step,  is  shown  in  Figure  4.10  to 
4.16.  Figure  4.10  is  a  plot  of  a  versus  y/h.  The  zero  flow  angle  is  indicated 
by  a  dashed  line  to  m^e  it  easier  to  identify  positive  and  negative  flow 
angles.  Figures  4.11,  4.12  and  4.13  are  of  the  non-dimensional  total  velocity 
and  non-dimensional  u  and  w  velocity  components  versus  y/h.  The  zero  velo¬ 
city  is  indicated  by  a  dashed  line.  Figures  4.14  and  4.15  are  the  two  velocity 
components,  u  and  w,  versus  the  ratio  of  probe  height  to  the  u  velocity  com¬ 
ponent  boundary-layer  thickness.  Again,  the  zero  velocity  is  indicated  by  a 
dashed  line.  The  final  plot.  Figure  4.16,  is  the  natural  log  plot  used  to  deter¬ 
mine  the  u  velocity  component  boimdary-layer  thickness.  The  constant  u 
component  edge  velocity  and  the  linear  curve  fit  is  shown  as  a  dashed  line. 
The  intersection  of  these  two  lines  indicate  the  boundary-layer  thickness. 

Measurement  Uncertainty 


Calibration  Data 

The  Magnehelic  pressure  gauges  used  to  calibrate  the  hot  film 
anemometer  h^  an  uncertainty  of  ±2%  full  scale,  according  to 
manufacturer’s  data.  At  the  lower  limit  of  each  gauge  this  imcertzunty 
translates  to  a  50%  error  in  determining  flow  velocity.  The  lowest  velocities 
read  from  each  gauge  were  as  follows:  for  the  0  to  10  inch  water  gauge,  68  fps; 
for  the  0  to  1  inch  water  gauge,  34  fps,  and  for  the  0  to  0.25  inch  water  gauge, 
9  fps.  It  was  felt,  based  on  the  extremely  od  lineeir  fit  of  the  calibration 
data  to  King’s  law,  that  the  actual  error  in  calculating  the  velocity  was  consid¬ 
erably  smaller.  Repeatability  of  the  calibration  data  was  within  ±0.5  fps 
which  corresponds  to  a  maximum  of  5%  difference  between  readings.  As 
stated  previously,  the  calibration  velocity  measured  from  the  pitot  tube  and 
the  velocity  calculated  from  the  fourth  degree  polynomial  yielded  a  maximum 
of  ±5%  difference.  From  information  described  in  TSI  Bulletin  TB5,  the  error 
due  to  a  fluid  temperature  difference  of  4  degrees  Fahrenheit  between  the  cali¬ 
bration  2md  velocity  profile  tests  is  considered  to  be  less  than  1%. 
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STEP  ANGLE-  15  degrees,  X-  1.00  Inches 
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Figure  4.10  Total  Flow  Angle  (a)  Versus  Figure  4.11  y/h  Versus  U/Ug  for 

y/h  for  A  =  15  Degrees  and  A  =  15  Degrees  and 

x/h  =  2.0  x/h  =  2.0 


STEP  RNGLE-  15  degrees,  X-  1.00  Inches 
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Figure  4.12  y/h  Versus  u/ug  for  Figure  4.13  y/h  Versus  w/Ug  for 

A  =  15  Degrees  and  A  =  15  Degrees  and 

x/h  =  2.0  x/h  =  2.0 
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Figure  4.14  y/6  Versus  u/Ug  for  Figure  4.15  y/6  Versus  w/Ug  for 

A=  15  Degrees  and  A  =  15  Degrees  and 


STEP  RNGLE-  15  degrees,  X-  1.00 


Ln(y) 


Probe  Positioning 


The  modified  height  gauge  of  the  probe  positioning  mechanism  had  a 
resolution  of  0.001  inches.  Because  of  the  problem  of  the  relative  movement 
of  the  model  to  the  wind  ttumel,  described  previously,  the  uncertainty  of  the 
probe  position  was  ±0.002.  The  repeatability  of  probe  positioning  was  very 
good  and  within  the  ±0.002  imcertainty.  The  probe  yaw  angle  was  set  by  an 
yaw  indicator  wheel  calibrated  from  0  to  360  degrees  in  1  degree  bicrements. 

Data  Reduction  -  Total  Velocity  and  Flow  Direction 

It  was  believed  that  the  data  reduction  procedure  could  resolve  the 
flow  direction  to  within  ±2.0  degrees.  From  repeated  velocity  profiles  at  the 
same  survey  locations,  the  flow  angle  showed  very  good  repeatability  when  the 
total  velocity  was  above  10  fps.  Below  10  fps  the  repeatability  of  the  flow 
angle  was  poor,  ±10  degrees,  but  this  was  to  be  expected  since  these  measure¬ 
ments  were  made  at  or  below  the  lower  limit  of  the  calibration.  The  uncer¬ 
tainty  of  the  data  reduction  procedure  in  determining  the  total  flow  velocity 
was  estimated  to  be  within  ±3  fps.  The  repeatability  of  the  calculated  totzd 
velocity  between  tests  at  the  same  survey  location  was  very  good  with  the 
velocities  differing  by  only  ±1  fps. 
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Chapter  5 

DISCUSSION  OF  RESULTS 
Overview 


The  final  part  of  this  study,  which  is  described  in  this  chapter,  was  to 
interpret  the  results  for  the  flow  over  a  swept  backward-facing  step. 
Observations  and  analysis  of  the  oil  flow  visualization  tests  are  discussed 
first.  Then  an  explanation  of  the  velocity  profile  survey  results  is  given.  An 
overview  of  the  flowfield  in  the  vicinity  of  each  step  is  presented  in  a  carpet 
plot  format.  A  more  detailed  look  at  the  recirculation  region  behind  each 
step  is  provided  2uid  discTissed.  A  comparison  of  the  test  results  from  this 
study  are  made  with  other  investigator’s  works.  The  results  outside  of  the 
separated  region  are  correlated  in  terms  of  Coles’  incompressible  law  of  the 
wall  and  wake.  Finally,  the  Baldwin-Lomax  algebraic  tmbulence  model, 
modified  for  swept  bacl^ard-facing  steps,  is  investigated.  An  attempt  is 
made  to  correlate  the  normal  (u)  and  parallel  (w)  velocity  components  in 
terms  of  parameters  of  the  Baldwin-Lomax  model. 

Oil  Flow  Results 


0  Degree  Step 

The  flow  visualization  studies  indicated  good  agreement  with  the 
results  presented  by  Selby  (1982).  Figure  5.1  is  a  typical  photograph  of  the 
oil  flow  behind  the  0  degree  step.  From  this  photograph  the  reattachment 
line  is  located  approximately  2.5  inches  or  5  step  heights  downstream  of  the 
step.  Selby’s  (1982)  results  indicate  reattachment  at  approximately  6.5  step 
heights  downstream  of  the  step.  The  actual  flow  directions  can  not  easily  ^ 
determined  from  the  photographs  but  by  observing  and  noting  the 
movement  of  the  oil  behind  the  step  during  the  tests  the  flow  directions 
could  be  resolved.  From  the  step  face  to  the  reattachment  line  the  oil  flowed 
back  towards  the  step,  and  accumulated  in  a  rather  large  puddle.  Because  of 
the  excessive  amoimt  of  oil  put  on  the  plate,  the  oil  accumulation  line 
extended  from  the  step  to  approximately  1.0  inches  downstream  of  the  step. 
It  is  believed  that  because  of  the  excessive  amoimt  of  oil  behind  the  step  the 
flowfield  was  disturbed,  but  the  extent  of  the  interference  is  unknown. 
Downstream  of  the  reattachment  line  the  oil  flowed  in  the  freestream 
direction.  The  flow  direction  inside  and  outside  of  the  recirculating  region 
appeared  to  be  perpendicular  to  the  face  of  the  step.  The  flow  visualization 
indicated  no  wall  interference  for  z  =±16.0  inches  about  the  centerline  of  the 
model.  The  reattachment  region  also  remained  parzdlel  to  the  step  fsure  over 
this  same  range. 

15  Degree  Step 

Photographs  showing  the  oil  flow  behind  the  15  degree  swept  step 
were  of  poor  quality  and  could  not  be  photocopied  to  show  any  flow  detail; 
therefore,  they  are  not  included.  From  the  observations  made  during  the 
testing,  the  reattachment  line  was  approximately  2.5  inches  (5.0  step 
heights)  in  the  streamwise  direction,  downstream  of  the  step  (2.4  inches  or 
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4.8  step  heights  measured  perpendicular  to  the  step).  Selby’s  (1982) 
reattachment  point  for  the  15  degree  step  was  approximately  x'/h=5.4 
downstream  of  the  step.  Between  the  step  face  and  reattachment  line,  small 
ripples  in  the  surface  of  the  oil  were  seen  moving  parallel  to  the  step  face. 
The  flow  in  this  region  was  also  moving  toward  the  step  as  evidenced  by  the 
accumulation  of  the  oil  near  the  step  face,  similar  to  the  result  obtained  in 
the  zero  degree  step  case.  Downstream  of  the  reattachment  line,  the  flow 
again  moved  in  the  direction  of  the  freestream.  Wall  effects  had  no  influence 
on  the  flow  between  z  =±12  inches  of  the  model  centerline.  The 
reattachment  line  remained  parallel  to  the  step  over  the  same  range  of  z. 

30  Degree  Step 

Figures  5.2  and  5.3  show  the  oil  flows  behind  the  30  degree  step. 

From  the  photographs  and  the  observations  the  reattachment  line  was 
determined  to  be  approximately  3.2  inches  (6.4  step  heights)  in  the 
streamwise  direction,  downstream  of  the  step  face,  or  2.8  inches  (5.6  step 
heights)  perpendicular  to  the  step.  Selby’s  (1982)  analysis  indicated  a 
reattachment  point  of  x'/h=5.4  downstream  of  the  step.  Figure  5.2  was 
taken  15  seconds  after  the  timnel  was  brought  up  to  the  test  velocity  of  50 
mph.  The  formation  of  ripples  moving  parallel  to  the  step  face,  and  the 
accumulation  of  oil  behind  the  step  can  be  clearly  seen.  Thus,  the  oil  was 
flowing  both  parallel  and  perpendicular  to  the  step  face.  Within  the 
recirculation  region  the  flc  v  was  trying  to  turn  so  that  it  was  perpendicular 
to  the  step.  Within  2.0  inches  downstream  of  the  reattachment,  the  flow  is 
s^ain  parallel  to  the  freestream  direction.  No  effect  due  to  wall  interference 
coxild  be  seen  in  the  flow  between  z  =±16  inches  of  the  centerline  of  the 
model.  The  reattachment  line  was  parallel  to  the  step  over  the  same  rmge 
of  z. 


Velocity  Profile  Results 


Overview 


The  test  results  for  the  three  backward-facing  step  configurations  are 
presented  and  discussed  in  the  following  sections.  The  presentation  of  these 
results  consists  of  two  sets  of  four  carpet  plots  each.  The  first  set  is  an 
overview  of  the  entire  range,  upstream  and  downstream  of  the  step,  with 
only  selected  profiles  being  shown  in  order  to  improve  the  clarity  of  the 
plots.  The  second  set  includes  all  of  the  profiles  within  the  recirculation 
region  of  each  step  configuration.  Each  set  was  made  up  of  the  following 
carpet  plot  profiles:  1)  total  velocity;  2)  total  velocity  flow  direction;  3)  u- 
component  of  the  total  velocity;  and  4)  w-component  of  the  total  velocity. 
To  further  improve  the  readability  of  the  carpet  plots,  the  individual  profile 
data  points  2U'e  connected  by  straight  lines  and  the  zero  position  of  each 
profile  is  indicated  by  a  dotted  line,  thus  making  it  easier  to  identify  positive 
and  negative  values. 
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Flow  Direction 
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the  30  Degree  Step,  Uod=  73  fps,  (15 
seconds  after  flow  began) 


Flow  Direction 


Degree  Step.Uoo"^  73  fps 
results) 


0  Degree  Swept  Step 

Overview  of  Velocity  Profiles.  Because  of  the  interest  in  the  flow  over 
a  backward-f£u:ing  unswept  step  an  abundance  of  experimental  results  are 
available.  Thus,  by  comparing  the  0  degree  step  results  with  previously 
published  data,  a  determination  of  the  validity  of  the  results  from  this  study 
can  be  made.  For  the  0  degree  step,  velocity  profiles  were  measured  at  the 
following  X  locations;  -2.75,  -1.5,  -0.5,  0.0,  0.1,  0.25,  1.0  to  4.0  in  steps  of  0.5, 
5.0,  6.0  and  7.5  inches.  All  survey  locations  were  measured  relative  to  the 
step,  with  negative  measurements  indicating  profile  locations  upstream  of  the 
step. 

Carpet  plots  of  the  0  degree  step  profiles  are  presented  in  Figures  5.4 
to  5.7.  Figure  5.4  shows  the  non-dimensional  total  velocity  profiles.  The 
total  velocity  flow  direction  (a)  carpet  plot  is  shown  in  Figure  5.5,  with  all 
flow  angles  being  approximately  0  degrees.  The  0  degree  flow  angle  implies 
that  the  flow  remains  in  the  free  stream  direction.  Figure  5.6  shows  tl  e  u 
velocity  components.  For  this  step  the  u  component  velocity  is  the  same  as 
the  total  velocity  because  of  the  near  0  degree  flow  angles.  The  w  velocity 
components  shown  in  Figure  5.7  are  all  practically  zero,  with  a  mayirmiTn 
value  of  w/u  e=0.08. 

Recirculation  Region  Velocity  Profiles.  The  area  of  most  interest  is 
the  recirculation  region  downstream  of  the  step.  Carpet  plots  of  all  velocity 
profiles  from  the  x/h=0.0  to  8.0  are  presented  in  Figure  5.8  and  5.9.  Only 
the  total  velocity  and  the  flow  angle  carpet  plots  are  shown.  The  u  and  w 
component  plots  provided  no  additional  information,  the  u  component  is  the 
same  as  the  total  velocity  and  the  w  components  were  all  nearly  zero.  The 
velocity  profiles  in  Figure  5.8,  at  x/h=:0.2  and  0.5,  decrease  to  nearly  zero 
velocity  just  below  y/h=1.0.  The  reattachment  of  the  flow  occurs  between 
x/h=4  and  5  and  because  of  the  large  negative  velocities  present  at  x/h=4 
the  reattachment  is  probably  closer  to  5.  This  result  is  in  good  agreement 
with  the  oil  flows  observations  of  this  study.  To  precisely  determine  the 
reattachment  point  more  STirveys  locations  would  be  necessary.  Downstream 
of  the  reattachment  region,  the  effect  of  the  step  on  the  profile  is  slowly 
dissipated.  At  the  reattachment  point  the  velocity  gradient  dU/dy  equals 
zero  but  as  can  be  seen  in  Figure  5.8,  dU/dy  increases  with  increasing 
downstream  distance.  For  the  most  part,  the  flow  angles.  Figure  5.9,  above 
y/h=1.0  are  within  the  range  of  ±4  degrees.  Below  y/h=1.0  zind  within  3 
step  heights  of  the  step,  the  flow  angle  has  a  range  of  approximately  ±20 
degrees.  The  results  for  the  two  profiles,  x/h=0.2  and  0.5  and  below 
y/h=1.0,  are  questionable  because  the  tot^  velocity  is  outside  of  the 
calibration  limits. 

Comparison  With  Other  Investigators’  Results.  A  comparison 
between  the  0  degree  step  results  of  this  study  auid  the  results  of  other 
investigators  was  made.  As  discussed  by  Eaton  and  Johnston  (1981),  they 
listed  five  principal  independent  parameters  that  afiected  the  reattachment 
point,  and  thus  the  flowfield  in  the  vicinity  of  a  backward-facing  step.  These 
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Figure  5.4  Carpet  PI 
for  A  =  1 
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Figure  5.5  Carpet  Plot  of  the  a  Profiles 

for  A  =  0  Degree,  x/h  =  -5.5  to  15 
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Figure  5.6  Carpet  Plot  of  the  u/u^  Profiles 
for  A  =  0  Degree,  x/h  =  -5.5  to  15 
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Figure  5.9  Carpet  Plot  of  the  a  Profiles 

In  the  Recirculation  Region  for 
A  =  0  Degree,  x/h  =  .2  to  8.0 


five  parameters  are:  l)  laminar  or  turbulent  preseparation  boundary-layer;  2) 
separation  boundary  layer  thickness;  3)  fireestream  turbulence;  4)  streamwise 
pressure  gradient;  and  5]  aspect  ratio  of  the  flow  apparatus  (channel 
width/step  height).  It  is  reasonable  to  assume  that  it  would  be  nearly 
impossible  for  two  investigators,  at  different  facilities,  to  match  all  five 
parameters;  consequently  slight  differences  in  the  flowfields  as  well  as  the 
restilts  must  be  expected.  The  amount  that  the  results  differ  will  depend  on 
how  close  the  parameters  are  matched. 

A  comparison  of  the  mean  velocity  profiles  of  Eaton  and  Johnston 
(1980),  Etheridge  and  Kemp  (1978)  and  this  study  are  presented  in  Figure 
5.10.  It  should  be  noted  that  the  Eaton  and  Johnston,  and  Etheridge  2nd 
Kemp  profiles  are  located  at  the  reattachment  point.  The  results  from  this 
study  represent  a  survey  location  near,  but  downstream  of  the  reattachment 
point.  It  also  must  be  noted  that  the  abscissa  of  the  other  investigators  are 
shifted  so  that  at  y/h=2  the  ratio  of  u/u  ^=1.  Takin;-,  ’uto  account  the 
differences  in  survey  locations,  the  agreement  between  the  results  is 
considered  good.  Bradshaw  and  Wong  (1972)  provided  a  mean  velocity 
profile  for  a  downstream  location  of  x/h=10.  The  mean  velocity  profile  from 
this  study  which  best  matches  these  results  is  located  at  x/h=15  and  the 
comparison  is  shown  in  Figure  5.11.  Once  again,  the  agreement  between 
these  results  is  very  good  with  a  percent  difference  of  less  than  5%.  From 
the  above  comparisons,  it  is  felt  that  because  of  the  good  agreement  with  the 
results  of  several  researchers  the  validity  and  accuracy  of  the  results  from 
this  study  are  acceptable,  especially  when  the  problems  involved  with 
matching  the  flowfields  are  taken  into  account. 

15  Degree  Swept  Step 

Selected  Velocity  Profiles.  There  are  no  known  experimental  results 
available  for  comparison  purposes  for  the  flow  over  a  swept  backward-fcicing 
step.  The  15  degree  step  velocity  profiles  were  measmed  at  the  following  x 
locations:  -2.56,  -1.5,  -0.5,  0.0,  0.15,  0.25,  0.5  to  4.0  in  steps  of  0.5,  5.0,  6.0 
and  7.5  inches.  Selected  profiles  of  the  15  degree  swept  step  are  presented  in 
Figures  5.12  to  5.15.  The  total  velocity  carpet  plot.  Figure  5.12,  indicates  a 
reattachment  point  between  x/h=4  and  6.  The  total  velocity  flow  angles. 
Figure  5.13,  are  nearly  zero  everywhere  except  within  and  just  downstream 
of  the  recirculation  region.  Figure  5.14  shows  the  u  velocity  component  is 
equal  to  the  total  velocity  except  in  and  just  downstream  of  the  recirculation 
region.  The  negative  u  velocity  component,  due  to  flow  reversal,  can  be  seen 
below  y/h=0.25  for  the  two  profiles  at  x/h=2  and  4.  The  w  velocity 
components.  Figure  5.15,  are  positive  everywhere  except  in  the  recirculation 
region. 


Recirculation  Region  Velocity  Profil^.  A  closer  look  at  the 
recirculation  region  reveals  more  details  of  the  flow  over  the  15  degree  swept 
backward-facing  step.  The  presence  of  the  step  greatly  affects  the  total 
velocity  profiles  at  x/h=0.3,  0.5  and  1.0,  Figure  5.16.  Just  below  the  step 
siirface,  y/h=1.0,  the  total  velocity  decreases  almost  to  zero.  The  profiles  at 


61 


15  DEG.  STEP,  U/Ue  VELOCITY  PROFILES 


63 


Figure  5.12  Carpet  Plot  of  the  U/Up  Profiles 
for  A  =  15  Degrees,  x/h  =  -5.1  t 
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Figure  5.13  Carpet  Plot  of  the  a  Profiles 
for  A  =  15  Degrees,  x/h  =  -5. 
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Figure  5.16  Carpet  Plot  of  the  U/Ug  Profiles  In  the  Recirculation  Region 
for  A  =  15  Degrees,  x/h  =  .3  to  8.0 


x/h=2.0,  3.0  and  4.0  indicated  a  strong  recirculation  region,  with  the 
stagnation  points  at  these  location  being  approximately  y/h=0.38,  0.3,  0.22, 
respectively.  Flow  reattachment  occurs  between  x/h=4  and  5,  but  probably 
closer  to  x/h=4.  More  survey  locations  are  required  to  accurately  determine 
the  reattachment  region.  Similar  to  the  0  degree  step,  the  velocity  gradient 
near  the  wall  increases  with  increasing  distance  downstream  of  the 
reattachment  point. 

As  is  expected  the  flow  angles.  Figure  5.17,  zire  nearly  zero  in  the 
freestream.  The  flow  angles  in  the  recirculation  region  and  near  the  surface 
downstream  of  the  reattachment  point  are  adl  negative.  For  a  given  profile 
location,  as  y/h  becomes  less  than  1.0,  the  flow  turns  so  that  it  is 
perpendicular  to  the  step  face.  The  flow  angles  for  the  three  profiles  at 
x/h=0.3,  0.5  and  1.0  and  y/h<1.0  are  questionable,  due  to  the  fact  that  the 
total  velocities  are  so  low  and  the  flow  reversal  ambiguity  could  not  be 
resolved.  The  next  three  profiles,  x/h=2.0,  3.0  and  4.0,  the  flow  reversal  was 
evident  and  the  profiles  have  been  modified  accordingly.  Once  below 
x/h=1.0,  the  influence  of  the  step  on  the  profiles  is  to  produce  a  negative 
flow  direction  angle.  One  important  point  that  should  be  noted  is  that 
although  the  flow  angle  is  negative  between  y/h  of  approximately  0.25  end 
the  surface  of  the  plate,  the  flow  is  actually  180  degrees  from  what  is  shown. 
This  is  due  to  the  flow  being  reversal  and  the  sign  of  the  total  velocity  being 
negative.  The  influence  of  the  step  is  still  present  in  the  flow  angle  data  at 
x/h=8.0,  but  the  effect  is  small  and  the  flow  is  gradually  realigning  itself 
with  the  freestream  flow  direction. 

The  u  velocity  components  are  shown  in  Figure  5.18.  The  u 
component  profiles  at  x/h=0.3,  0.5  and  1.0  decreases  rapidly  jiist  below 
y/h=1.0  because  of  the  presence  of  the  step.  Profiles  at  x/h=2.0,  3.0  zind  4.0 
indicates  a  strong  recirculation  region,  while  profiles  beyond  the 
reattachment  point  continue  to  show  some  influence  due  to  the  step. 

The  w  velocity  component,  Figure  5.19,  at  x/h=0.3,  0.5,  zmd  1.0  axe 
essentially  zero  below  y/h=1.0.  Profiles  between  and  including  x/h=2.0  and 
7.0  show  the  w  component  decreasing  from  the  free  stream  value  to  zero  or 
slightly  negative  and  then  increase  again  before  becoming  zero  at  the  surface. 
The  positive  w  velocity  component,  indicated  in  the  recirculation  region  near 
the  surface  of  the  plate,  is  responsible  for  producing  the  ripples  moving 
parallel  to  the  step  face  which  were  observed  in  the  oil  flow  studies.  Beyond 
x/h=7,  the  influence  of  the  step  is  small,  with  the  w  component  profiles 
remaining  constant  within  the  free  stream  region  but  then  decreasing  to  zero 
below  y/h=1.0.  The  reason  for  the  decrease  in  the  w  component  with 
decreasing  y/h  is  because  the  total  velocity  is  also  decreasing  to  zero.  ^ 
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Figure  5.18  Carpet  Plot  of  the  u/ug  Profiles  In  the  Recirculation  Region  for 
A  =  15  Degrees,  x/h  =  .3  to  8.0 
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Region  for  A=  15  Degrees,  x/h 


30  Degree  Swept  Step 

Selected  Velocity  Profiles.  The  30  degree  step  was  the  most 
thoroughly  tested  configuration  of  the  three.  The  same  number  of 
downstream  locations  was  surveyed,  but  a  more  complete  survey  upstream  of 
the  step  was  performed  in  order  to  more  accurately  define  the  initial  flow 
conditions.  The  following  profile  locations  were  srjveyed:  -2.5  to  0.0  in  steps 
of  0.5,  0.15,  0.25;  0.5  to  4.0  in  steps  of  0.5,  5.0,  6.0  and  7.5  inches. 

The  overview  of  the  total  velocity  profiles.  Figure  5.20,  reveals  very 
little  information.  The  upstream  profiles  are  typical  of  the  other  two 
conngurations.  Downstream  of  the  step,  in  Figure  5.20,  no  discernible 
recircTilation  region  can  be  observed,  but  when  all  the  profiles  downstream  of 
the  step  are  examined  the  recirctilation  region  becomes  apparent.  The  flow 
reattachment  point  is  believed  to  occur  between  x/h=2  and  4.  Downstream 
of  the  reattachment  point,  the  flow  slowly  dissipates  the  effect  of  the  step, 
with  the  influence  of  the  step  still  observable  at  x/h=:15.  The  totaJ  flow 
angle,  Figure  5.21,  is  basically  zero  at  the  following  profile  locations;  x/h=-5, 
-3,  0,  12  and  15.  At  the  remaining  survey  locations,  the  flow  angle  is  zero  in 
the  freestream  region  but  below  y/h=1.0  the  flow  angle  becomes  negative. 
The  overall  u  velocity  component  data,  Figure  5.22,  are  similar  to  the  total 
velocity  profiles.  Figure  5.23  displays  the  w  velocity  component  carpet  plots. 
The  w  components  are  positive  everywhere  except  below  y/h=1.0  and 
between  x/h=2  and  6  where  the  component  become  slightly  negative.  A 
better  view  of  the  flow  in  the  vicinity  of  the  30  degree  step  can  be  seen  by 
examining  the  profiles  between  the  step  and  x/h=8. 

Recirculation  Region  Velocity  Profiles.  The  total  velocity  profiles, 
Figure  5.24,  shows  a  recirculation  region  with  the  reattachment  region 
between  x/h=3.0  and  4.0.  More  survey  positions  are  necessary  to  accurately 
locate  the  reattachment  region.  The  profile  at  x/h=2  has  not  been  reversed 
because  a  definite  flow  reversal  could  not  be  detected  from  the  total  velocity 
profile.  It  is  felt  that  the  flow  did  reverse  at  this  location,  but  it  is  not 
discernible  from  the  measurements  .  The  profiles  at  survey  locations, 
x/h=0.3,  0.5  and  1.0  are  greatly  affected  by  the  close  proximity  of  the  step. 
One  noticeable  difference  is  that  the  decrease  in  totsd  velocity  is  not  as  great 
as  for  the  same  profiles  of  the  0  and  15  degree  configurations.  The  total 
velocity  ratio  below  y/h=l  for  the  0  and  15  degree  step  was  approximately 
0.07,  where  as  the  ratio  for  the  30  degree  step  was  0.20  almost  three  times 
larger.  As  stated  before,  the  results  of  this  region  are  questionable  due  to 
the  limits  of  the  calibration  and  the  lack  of  the  ability  to  resolve  the  flow 
reversal  ambiguity.  Downstream  of  the  reattachment  region,  the  surface 
velocity  gradient,  dU/dy,  increases  with  increasing  strezunwise  distance. 

The  total  flow  angles,  Figure  5.25,  are  approximately  zero  in  the 
freestream,  but  become  largely  negative  in  the  recirculation  region.  The 
three  survey  locations  of  x/h=0.3,  0.5  and  1.0  below  y/h=1.0  have  a 
constant  flow  angle  of  approximately  -50  degrees.  It  must  be  noted,  that 
because  of  the  choice  of  fl=15  degrees,  if  |a|>45  degrees  then  this  particular 
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Figure  5.20  Carpet  Plot  of  the  U/Ue  Profiles 
for  A  =  30  Degrees,  x/h  =  -5  to  15 
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Figure  5.21  Carpet  Plot  of  the  a  Profiles 

for  A  =  30  Degrees,  x/h  =  -5  to  15 
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Figure  5.24  Carpet  Plot  of  the  U/«e  Profiles  in  the  Recirculation 
Region  for  A  =  30  Degrees,  x/h  =  .3  to  8.0 


Figure  5.25  Carpet  Plot  of  the  a  Profiles  In  the  Recirculation 
Region  for  A  =  30  Degrees,  x/h  =  .3  to  8.0 


a  is  only  an  estimate.  This  is  a  requirement  of  the  yaw  calibration  curve, 
which  was  discxissed  in  the  ''’’aw  Calibration  section  of  Chapter  4.  At  survey 
locations,  x/h=2,  3  and  4,  and  below  y/h=1.0,  the  flow  angle  becomes 
negative  decreasing  to  a  minimxim,  then  increasing  towards  zero  degrees  and 
finally  decreasing  again.  For  x/h=2,  a  minimum  flow  angle  of  -50  degrees  at 
y/h=0.5  is  achieved.  As  y/h  continues  to  decrease  the  flow  angle  increases 
to  a  maximum  within  the  shear  layer  of  -5.3  degrees  at  y/h=0.08  and  then 
decreases  to  -9.3  at  y/h=0.04.  The  survey  location  at  x/h=3,  yields  a 
minirmiTn  flow  angle  of  -44  degrees  at  y/h=0.4  which  increases  to  a 
maximum  of  -4  degrees  at  y/h=0.12  and  then  finally  decreases  to  -26  at 
y/h=0.04.  The  stagnation  point  for  each  profile  occurs  at  the  maximum  flow 
angle  within  the  shear  layer.  Below  the  stagnation  point,  the  flow  angle 
again  decreases,  but  now  the  total  velocity  is  negative,  meaning  the  flow  is 
moving  toward  the  step.  For  x/h=4  the  mminuim  flow  angle  of  -38  degrees 
occurs  at  y/h=0.40  and  then  increases  to  -i-7  degrees  at  y/h=0.04.  The 
remaining  locations,  x/h=6,  7  and  8  continues  to  show  a  strong  influence  on 
the  flow  angle  due  to  the  presense  of  the  step.  By  location  x/h=8  the  effect 
of  the  step  is  beginning  to  dissipate.  For  profiles  x/h=6,  7  and  8  the 
minimum  flow  angle  occurs  much  closer  to  the  stirface  of  the  plate  than  the 
profiles  further  upstream. 

The  u  velocity  component  profiles.  Figure  5.26,  are  similar  to  the 
total  velocity  profiles.  One  noticeable  difference  is  that  the  transition  from 
the  freestream  to  the  shear  layer  is  not  as  smooth  for  the  u  component  as  for 
the  total  velocity. 

The  more  interesting  results  are  the  w  velocity  component.  Figure 
5.27.  The  carpet  plot  indicates  a  positive  w  component  outside  of  the  shear 
layer  with  the  component  becoming  negative  below  y/h=l  at  survey 
locations  x/h=0.3,  0.5  and  1.0.  The  profiles  at  x/h=2,  3  zuid  4  and  within 
the  shear  layer  show  the  w  velocity  component  suddenly  decreasing  and 
becoming  negative.  The  minimum  velocity  for  the  three  profiles  occurred  at 
y/h=0.5,  0.4  and  0.4,  respectively.  The  w  component  then  increases  to  a 
positive  value  as  the  surface  of  the  plate  is  approached.  This  supports  the 
observations  made  during  the  oil  flow  studies  which  indicated  that  the  oil 
flowed  in  the  positive  w  direction  within  the  recirculation  region  (see  Figure 
5.2).  Once  the  w  components  at  locations  x/h=5,  6  and  7  enter  the  shear 
layer  they  become  negative  and  remain  negative  to  the  surface  of  the  model. 
For  the  x/h=8  the  w  velocity  components  are  positive  along  the  entire 
profile. 

Comparison  of  Results 

Initially  there  was  a  concern  that  because  of  the  three  different 
leading  edge  eind  step  configurations,  identical  flowfields  could  not  be 
reproduced.  However,  the  reduced  data  show  that  the  total  velocity  profiles 
at  x/h=0.0  agree  to  within  about  4%  for  the  three  configurations.  From 
Schlichting  (1979)  the  following  relation  V3s  used  to  estimate  the  turbulent 
boundary-layer  thickness  at  the  edge  of  the  step: 
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Figure  5.27  Carpet  Plot  of  the  m/u^  Profiles  In  the  Recirculation 
Region  for  A  =  30  Degrees,  x/h  =  .3  to  8.0 
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where  u  =  50  mph  or  73.3  fps,  u  =  1.57  x  10“^  ft^/sec  and  Ljep  =  2.25  ft. 

This  equation  yields  a  £  of  0.052  ft  or  0.62  inches.  The  agreement  between  the 
experiment2Ll  and  theoretical  is  good  with  an  approximate  5%  difference. 

Using  the  procedure  described  in  Chapter  4,  the  boundary-layer  thick¬ 
ness  was  estimated  for  each  profile  surveyed.  A  plot  of  6  verses  x'  for  each 
swept  step  configuration  is  shown  in  Figure  5.28.  The  results  upstream  of  the 
step  indicate  a  growing  boundary-layer  thickness  as  the  step  is  approached, 
with  6  +  h  =1.15  inches  at  x/h=0.0.  It  should  be  noted  that  the  boundary- 
layer  thickness  is  indicated  with  respect  to  the  lower  surface.  Thus,  the  S ’s 
upstream  of  the  step  are  0.5  inches  thicker  to  accoimt  for  the  height  of  the 
step.  The  furthest  upstream  result  for  the  0  degree  step  is  questionable. 
Downstream  of  the  step,  x/h=0.0  to  6.0,  S  decreases  to  approximately  0.70 
inches.  Beyond  the  recirculation  region,  x/h  >  6.0  the  boundary-layer  thick¬ 
ness  again  begins  to  increase  as  the  flow  recovers  from  the  separation  caused 
by  the  step.  In  the  recovery  region,  S  generally  increases  with  increasing 
sweep  angle  for  a  specific  survey  position. 

Estimation  of  Stagnation  Line  Within  the  Recirculation  Region.  An 
attempt  was  made  to  determine  the  position  of  the  stagnation  line  for  both 
the  u  and  w  velocity  components  within  the  recirculation  region.  For  the  u 
component,  this  was  accomplished  by  examining  the  total  velocity  profiles 
that  showed  evidence  of  flow  reversal.  Then,  by  linear  interpolating  between 
the  two  data  points  on  either  side  of  the  u  component  stagnation  line,  a  value 
of  y  when  u=0  was  determined.  The  results  are  presented  in  Figure  5.29.  The 
results  for  the  0,  15  and  30  degree  step  au'e  fitted  to  separate  second  degree 
polvnomials  using  a  least  squares  method.  The  results  for  the  0  degree  step  at 
X/  1=2.0  amd  the  15  degree  step  at  x/h=3.0  and  4.0  are  neglected  because  it  is 
not  clear  where  the  flow  reversal  occurs.  However  from  the  curve  fits,  an  esti¬ 
mate  of  the  location  of  the  reverse  flow  region  can  be  made.  The  general 
trend  is  that  the  reattachment  length  decreases  as  the  sweep  angle  increases. 

The  w  component  stagnation  line  was  found  by  noting,  for  each  survey, 
the  height  above  the  wall  where  w=0  for  the  first  time  as  y/h  decreased  from 
4.0  to  0.0.  Since  the  w  component  is  essentially  zero  for  the  0  degree  step, 
only  results  from  the  15  and  30  degree  steps  are  shown  in  Figure  5.30.  From 
these  results  the  stagnation  line  for  w  component  is  linear.  Figure  5.30  indi¬ 
cates  that  within  the  reattachment  region  (where  u=0)  the  w  component  still 
has  finite  magnitude.  The  w  component  is  also  negative  at  this  location  which 
means  the  external  flow  has  been  turned  beyond  normal  to  the  step  by  the 
pressure  gradient  acting  normal  to  the  step,  so  that  at  reattachment  the  flow 
is  moving  in  the  negative  z'  direction.  This  is  consistent  with  the  oil  flow 
observation,  presented  in  Figure  5.3. 
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Boundary— Layer  Thicknesses 
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Figure  5.28  Boundary-Layer  Thickness  as  a  Function 
of  x'/h  for  the  Three  Step  Sweep  Angles 


Downstream  Height 
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Line  Within  the  ItecTrculatlon  Region  for  the 
Three  Swept  Step  Configurations 


Correlation  of  Results  to  Computational  Turbulence  Models 
Overview. 


Turbulent  boundary-layers  are  made  up  of  two  regions:  zui  inner  or  wall 
region  and  a  outer  or  wake  region.  The  inner  region,  which  consists  of  only 
alwut  10  to  20%  of  the  total  turbulent  boundary-layer,  may  be  further  broken 
down  into  three  sublayers.  The  layer  closest  to  the  wsill  is  the  laminar  or 
viscous  sublayer  where  the  viscous  shear  dominates.  This  layer  extends  from 
the  wall  to  approximately  y"*"  =  5,  where  y"*"  =  y-a^ju.  The  buffer  sublayer, 
where  viscous  and  turbulent  stresses  are  approximately  equal,  extends  from 
5<  y*^  <35.  The  turbulent  stresses  are  dominant  in  the  final  sublayer  of  the 
inner  region,  the  fully  tiirbulent  layer,  which  extends  from  35<y‘'‘  to  y/6«0.2. 
The  outer  or  wake  region  makes  up  approximately  80%  to  90%  of  the  tur¬ 
bulent  botmdary-layer.  This  layer  is  characterized  by  the  intermittency  of 
turbulent  flow.  The  outer  region  is  more  sensitive  to  external  parameters, 
especially  the  pressure  gradient.  The  main  reason  for  needing  an  inner  and 
outer  region  to  numerically  approximate  turbulent  boundary-layers  is  because 
they  are  based  on  different  length  scales.  The  inner  region  length  sczde,  ujur, 
depends  on  the  wall  shear  stress  and  the  outer  region  length  scale  is  the  total 
boundary-layer  thickness,  8. 

One  aspect  of  this  study  was  to  investigate  the  correlation  of  the 
results  with  the  widely  used  turbulent  boundary-layer  law  of  the  wall  and  law 
of  the  wake  to  swept  backward-facing  step  flows. 

Coles’  Law 


Description  of  Law.  The  first  approximation  examined  was 
Coles’  incompressible  law  of  the  wall  and  wake  for  attached  flows  as  described 
by  Cebeci  and  Smith  (1974).  Coles’  equation  is  written  as: 


u+  =  -ln(y+)  +  C  + 


2n  .  2 

- sm 

K 


UL 

28 


(5.2) 


where  u"*"  =  u/u  ,  and  u,.  =  Vr/p  .  Coles’  law  is  valid  for  the  fully  turbulent 
sublayer  and  across  most  of  the  outer  region.  The  first  two  terms  on  the  right 
hand  side  of  the  equation  represent  the  semi-logarithmic  form  of  the  law  of 
the  wsJl.  K,  the  von  Karman  constant,  is  the  inverse  of  the  slope  and  C  is  the 
intercept  of  the  semi-logzu-ithmic  law  of  the  wedl  curve.  The  last  term  is 
Coles’  wake  function,  which  accounts  for  the  wake  defect  of  the  velocity  in  the 
outer  region  of  the  boundary-layer. 

The  four  constants  in  (5.2)  are  determined  for  the  present  data  by  a 
least  squares  technique.  This  is  accomplished  by  rewriting  Coles’  equation  in 
the  following  form  which  was  used  by  Danberg  (1971); 

u/u  e  =  ai  ln(y)  +  aj  +  as  sin^  ( a4  y )  (5.3) 
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The  a’s  are  four  profile  constants  to  be  determined  by  the  least  squzures  fit.  A 
general,  nonlinear  least  squares  fit  routine  was  used  to  compute  the  four  con¬ 
stants  from  four  initial  estimates  by  assuming  small  perturbations  to  the  a’s. 
This  allowed  the  equations  to  be  expressed  by  a  linear  Taylor  series  expansion, 
whereby  an  improved  estimate  of  the  a’s  was  obtained  to  start  a  new  itera¬ 
tion.  The  procedure  was  repeated  until  a  predetermined  degree  of  accuracy 
was  achieved.  With  the  four  values  of  a’s  known,  the  constants  Ur/u  g,  C,  11, 
and  6  could  be  determined.  It  is  not  possible  to  determine  k  independently 
and  the  normal  value  of  k  =  0.43,  used  by  Danberg  (1971),  was  assiuned. 

Results  from  Coles’  Law.  After  applying  the  procedure  described  above 
to  the  results  of  this  study,  up  and  downstream  of  the  recirculation  region,  the 
foiir  constants  u,./u  g,  C,  11  and  S  were  determined.  Instead  of  Ur/u  e  skin 
friction  coefficient,  Q,  was  calculated.  Cf  is  related  to  u^/u  g  by  the  following 
relation: 

'iS 

Cf  =  2  —  (5.4) 

V  / 

The  values  of  Cf  for  the  three  conffgurations  zu'e  plotted  in  Figure  5.31.  The 
Cf’s  upstream  of  the  step  indicate  the  presence  of  a  pressure  gradient  acting 
on  the  three  different  step  configurations.  For  the  0  degree  step  the  Cf’s 
decrease  as  they  approach  the  edge  of  the  step  indicating  an  adverse  pressure 
gradient  acting  on  the  flow.  A  favorable  pressure  gradient  acts  upstream  of 
the  15  degree  step  with  Cf  increasing  as  the  step  is  approached.  The  30 
degree  case  indicates  a  zero  pressure  gradient  or  a  slight  favorable  pressure 
gradient.  Cf’s  upstream  of  the  30  degree  step  remain  nezu'ly  constant.  Down¬ 
stream  of  reattachment  the  Cf’s  for  the  0  degree  step  are  approximately  a  fac¬ 
tor  of  3  smaller  than  the  values  obtained  upstream.  The  Cf’s  for  the  15  and  30 
degree  steps  are  only  slightly  lower  downstream  when  compared  to  the  vadues 
upstream.  At  reattachment  the  values  of  Cf  is  0  because  the  velocity  gradient 
at  the  wall  is  0.  As  x'/h  increases  from  reattachment  to  approximately  5  step 
heights  downstream,  Cf  adso  increases  and  then  remains  nearly  constant  with 
further  increases  of  x'/h.  The  effect  of  sweep  angle  on  the  vaJue  of  Cf  indi¬ 
cates  that  as  sweep  angle  increases  so  does  Cf. 

Figure  5.32  shows  the  various  values  of  the  intercept  C  as  a  function  of 
survey  location  for  various  sweep  angles.  Upstream,  C  for  the  0  degree  step, 
under  an  adverse  pressure  gradient,  increases  with  increasing  x'/h.  The 
upstream  C’s,  for  15  degree  step,  decrease  ^  the  step  is  approached,  indicat¬ 
ing  a  favorable  pressure  gradient.  For  the  30  degree  step  with  a  zero  or 
slightly  favorable  pressme  gradient  the  values  of  C  remain  nearly  constant. 
Downstream  of  reattachment  the  C’s  tend  to  decrease  with  increasing  sweep 
angle.  The  value  of  C  is  zero  at  reattachment  and  then  increases  with  increas¬ 
ing  downstream  distance.  A  typical  value  of  C  according  to  Cebeci  auid  Smith 
(1974)  is  5.0  for  a  two-dimensional  flat  plate.  The  compeirison  between  their 
value  and  the  value  of  C  from  this  study  shows  the  effect  of  separation  and 
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SKIN  FRICTION  COEFFICIENT 


INTERCEPT, 


INTERCEPT  IN  LAW  OF  THE  WALL,  C 


x'/h 


Figure  5.32  The  Law  of  the  Wall  Intercept,  C,  as  a 
Function  of  x'/h  for  the  Three  Step 
Sweep  Angles 


89 


the  three-dimensional  flow  caused  by  sweep. 


The  wake  parameter,  IT,  is  shown  in  Figure  5.33,  for  the  three  sweep 
angles  and  the  different  survey  locations.  11  is  a  function  of  the  pressure  gra¬ 
dient  acting  on  the  viscous  flowfield.  Coles  (1968)  evaluated  a  large  ntimber 
of  incompressible  turbulent  boimdary-layer  profiles  on  flat  plates  at  high  Rey¬ 
nolds  numbers  (Re  >  3000)  to  obtain  a  value  for  11  of  0.55.  In  the  present 
case,  upstream  of  the  step,  an  average  value  of  0.71  for  IT  was  determined.  If 
the  two  highest  values  are  neglected,  the  average  11  is  0.61.  Upstream  of  the 
step,  the  trends  for  the  various  sweep  angles  can  be  better  seen  in  the 
expanded  scale  of  Figure  5.34.  The  general  trend  of  IT  is  very  similar  to  the 
trend  of  C  in  Figure  5.32  and  is  consistent  with  the  skin  friction  coefficient 
data,  n  increases  with  increasing  x  /h.  for  the  0  degree  step.  For  the  15 
degree  step,  11  decreases  as  the  step  is  approached,  and  in  the  30  degree  case 
n  remains  nearly  constant  or  decreases  slightly. 


Downstream  of  reattachment  (Figure  5.33)  the  wake  is  the  dominant 
component  of  the  boundary-layer.  The  inner  or  wall  layer  effectively  disap¬ 
pears  at  reattachment.  Coles’  profile  equation  becomes, 


u 


2u,.n  . 

- sm 

u^/c 
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(5.5) 


where  in  the  limit  as  tends  towards  zero  at  reattachment  the  term 


2u,.n 

lun - 

r,-*0  Ue  K. 


=  1 


(5.6) 


Thus  n  goes  to  infmity  as  u^  goes  to  zero.  This  is  only  true  in  the  limit  and 
as  the  boundary-layer  develops  downstream  of  reattachment  the  effects  of  the 
wall  layer  become  important  again.  This  result  is  evident  in  the  calculated 
values  of  11  downstream  of  reattachment.  Near  reattachment  11  is  quite  Izu-ge 
with  values  eiround  10.0.  As  the  downstream  distance  increases  the  magnitude 
of  n  decreases' to  3.45,  1.89  and  1.01  for  0, 15  and  30  degree  sweep  angles, 
respectively.  Although  11  is  slowly  decreasing  it  is  still  much  larger  than  its’ 
pre-step  vjilue.  The  results  also  indicate  that  IT  decreases  as  sweep  angle 
increases. 

The  fmad  parameter.  Coles’  6,  is  plotted  in  Figure  5.35.  Also  shown  in 
this  figure  are  the  boimdary-layer  thicknesses  from  Figure  5.28.  The  agree¬ 
ment  between  the  two  procedures  for  determining  the  boimdary-layer 
thicknesses  is  quite  good  with  the  general  trend  of  the  results  being  similar 
between  the  two  procedures.  The  general  trend  of  Coles’  6,  upstream  of  the 
step  for  all  sweep  angles,  is  to  increase  as  the  step  is  approached,  and  except 
for  a  few  points,  the  agreement  is  good.  No  5’s  are  determined  from  Coles’ 
law  in  the  recirculation  region.  In  the  recovery  region  the  boimdary-layer  con¬ 
tinues  to  grow  with  increasing  x'/h.  The  effect  of  sweep  angle  on  S  in  the 
recovery  region,  is  that  for  a  specific  location,  6  increases  with  increasing  A. 
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LAW  OF  THE  WAKE  PARAMETER,  U 


x'/h 


Figure  5.33  Coles'  Wake  Parameter  as  a  Function 
of  x'/h  for  the  Three  Step  Sweep 
Angles 
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Boundary— Layer  Thicknesses 
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Figure  5.35  Comparison  of  Coles'  Law  Boundary- Layer  Thickness 
and  Those  Determined  Directly  From  the 
Experimental  Results 


Another  way  of  examining  the  results  from  the  Coles’  law  analysis  is  by 
looking  at  the  semi-log  profiles.  The  semi-log  profiles  at  the  step  (x'/h  =  0.0) 
are  shown  in  Figure  5.36  for  the  three  A’s.  As  pointed  out  earlier,  the  profiles 
u/u^  versus  ytS  show  very  little  differences,  «  4%,  but  the  log  profiles  show 
considerable  aifference  because  they  present  the  profiles  in  a  more  sensitive 
manner.  The  relative  magnitude  of  the  wake  defect  is  consistent  with  the 
previous  discussion  of  the  coefficients  of  Coles’  law.  The  A  =  0  degree  profile 
is  characteristic  of  an  adverse  pressure  gradient  because  of  the  large  wake 
defect,  whereas  the  A  =  15  degrees  profile  shows  evidence  of  a  favorable 
pressure  gradient.  The  A  =  30  degrees  also  has  a  slight  favorable  pressure 
gradient  but  is  more  representative  of  an  equilibrium  flow.  It  is  not  clear 
whether  these  gradients  are  entirely  due  to  initial  conditions  or  if  the  sweep 
induces  them.  Figure  5.37  presents  three  typical  semi-log  profiles  zdter  flow 
reattachment  for  A  =  30  degrees.  The  three  profile  locations  are  x/h=6.0,  7.0 
and  10.0.  These  three  profiles  illustrate  the  recovery  of  the  flow  after  reat¬ 
tachment.  The  wake  effect  is  quite  large  in  the  mofile  at  x/h=6.0.  At 
x/h=7.0  the  wake  effect  has  been  reduced  by  50%  and  at  x/h=10.0  the  wake 
effect  is  approximately  25%  of  what  it  was  at  x/h=6.0. 

The  resulting  velocity  profiles  produced  from  Coles’  equations  are 
shown  in  Figure  5.38  through  5.41  along  with  the  experimental  results  of  this 
study.  Since  the  upstream  velocity  profiles  are  similar  for  the  three  configura¬ 
tions  only  the  A  =  15  degrees  results  are  presented  in  Figure  5.38.  All  profiles 
upstream  of  the  step  collapse  into  approximately  a  single  curve,  with  the  solid 
lines  representing  Coles’  law  and  the  symbols  indicating  the  experimental 
data.  Figures  5.39  to  5.41  show  the  profile  from  Coles’  law  downstream  of 
reattachment.  All  A’s  show  recovery  from  separation.  At  separation  the  loga¬ 
rithmic  terms  of  Coles’  equation  (5.3)  are  zero  and  the  separation  profile  is 
just  the  wedce  component.  As  the  downstream  distzmce  increases  the  profiles 
become  fuller  indicating  an  increasing  importance  of  the  logau-ithmic  terms. 

Baldwin-Lomax  Method 


General  Baldwin-Lomax  Method.  The  method  as  described  by  Baldwin 
and  Lomax  (1978)  is  a  two-layer  algebraic  turbulence  model  which  is  widely 
used  in  computational  fluid  dynamics.  The  model  works  well  for  two-  and 
three-dimensional  attached  flows.  Basically  the  model  evaluates  the  turbulent 
shear  stresses  in  terms  of  an  eddy  viscosity.  The  major  advantage  of  this 
model  is  that  the  boundary-layer  thickness  does  not  have  to  be  known  expli¬ 
citly  and  thxis  the  arbitrariness  and  potential  error  associated  with  this  quanti¬ 
ty  is  eliminated. 

In  the  Baldwin-Lomax  model  the  eddy  viscosity,  /it,  is  represented  by 
two  values,  one  for  the  inner  and  one  for  the  outer  regions.  The  iimer  region 
uses  the  Prandtl-Van  Driest  formulation  for  /if  For  this  study,  only  the  outer 
region  formulation  was  used,  since  no  measurements  were  obtained  within  the 
inner  region.  The  value  of  the  outer  region  eddy  viscosity  coefficient  depends 
on  whether  the  flow  is  attached  or  separated.  For  attached  flow  the  outer 
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LAW  OF  THE  WALL  -  WAKE  PROFILE 

x'/h  =  0.0 


LOG(Y  +  ) 


Figure  5.36  Law  of  the  Wall  Semi -Logarithmic  Wake 
Profiles  at  x'/h  *  0.0  for  the  Three 
Step  Sweep  Angles  (symbols  are 
experimental  data) 
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LAW  OF  THE  WALL  -  WAKE  PROFILE 

AFTER  REATTACHMENT 
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Figure  5.37  Law  of  the  Hall  Serai -Logarithmic  Wake  Profiles 
After  Reattachment  for  A  =  30  Degrees 
(syndrals  are  experimental  data) 
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Figure  5.41  Coles’  Law  Velocity  Profile  After 
Reattachment  for  A  =  30  Degrees 
(symbols  are  experimental  data) 


eddy  viscosity  is  given  by; 

A*t  ~  KCcpPyjQax  ^max  ^  (5.7) 


where  K  is  the  Clauser  constant  equal  to  0.0168,  C^p  is  a  constant  equal  to 
1.6,  p  is  the  density,  ymax  is  the  nonnad  distance  to  Fmax  is  deter¬ 

mined  from  the  moment  of  vorticity  equation: 


F(y)=y|(*;|  1  —  e^ 


(5.8) 


where  u;  is  the  vorticity  and  A"*"  is  a  constant  equal  to  26.  ^  is  the 
KlebanofTs  intermittency  factor  defined  by  the  following  equation: 


Ckleby 

6‘ 

1  +  5.5 

ymax 

V  > 

(5.9) 


where  Cygb  = 


The  outer  layer  eddy  viscosity  for  a  separated  flow  is  given  by  the  fol¬ 
lowing  relation: 

Ud  Ymax 


=  K  p  Cep 


max 


(5,10) 


Ud  =  vTu^  +  v7  +  w7)|^  - 


(5.11) 


where  the  constant  =  0.25 


is  the  difference  between  the  maximum  and  minimum  velocity  across  the 
profile  and  Ujpin  =  0  for  attached  flow  due  to  the  no  slip  boundary  condition. 


Qualification  to  the  Baldwin-Lomax  Model.  In  order  to  achieve  satis¬ 
factory  resuItsTwinrnie'BalHwnrLomajrturSuIenw  model  several  qualifica¬ 
tion,  presented  by  Danberg  and  Patel  (1988),  need  to  be  made  to  the  model. 
First,  the  equation  defining  the  eddy  viscosity  for  the  separated  or  wake  flow 
case  must  be  modified.  The  necessary  changes  to  equation  (5.8),  which  de¬ 
fines  the  moment  of  vorticity,  F(y),  is  given  as  follows: 


F(y)  =  (y  -  hy)  (du/dy) 


1 


g(-y+/A  +) 


(5.12) 


The  first  difference  between  (5.12)  and  (5.8)  is  that  ui  (5.12),  the  quan¬ 
tity  (y  —  hy)  has  replaced  y,  with  hy  being  definea  as  the  vertical  height 
above  the  surface  to  the  position  of  the  minimum  velocity,  see  Figure  5.42. 
The  purpose  of  (5.8)  or  (5.12)  is  to  determine  the  appropriate  length  and 
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NOTE:  For  typical  attached  profiles  yg,ax  »  .65  5 

Figure  5.42  Typical  Separated  Profile  Illustrating  Terms  Important 
to  Modification  to  the  Baldwin- Lomax  Model 
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velocity  scale  for  evaluating  the  eddy  viscosity.  If  the  equation  is  left  un¬ 
changed,  as  stated  in  (5.8),  then  an  incorrect  value  of  y^ax  is  obtained.  This 
can  be  seen  by  considering  two  typical  profiles;  one  just  upstream  of  separa¬ 
tion  and  the  other  immediately  ^ter  separation.  In  the  first  case,  y^j^x  is 
roughly  65%  of  the  boimdary-layer  thickness,  which  is  characteristic  of  the 
thickness  of  the  wake  layer.  Note  that  the  maximum  vorticity  occurs  at  the 
wall  but  the  moment  arm  is  zero  there,  and  thus  the  maximum  of  the  mo¬ 
ment  of  vorticity  occurs  away  from  the  wall  at  a  point  characteristic  of  the 
wake  region. 

Immediately  after  separation  the  maximum  vorticity  still  occurs  near 
the  top  of  the  step,  but  now  the  corresponding  moment  arm  is  large  and  the 
maximum  in  F(y)  occurs  at  that  point.  The  characteristic  length  scale  sud¬ 
denly  Jumps  from  65%  of  S  to  the  step  height,  H.  As  a  result  there  is  a 
discontinuity  in  the  length  scale  y^ax  Fmax»  which  is  physically  unrealis¬ 
tic.  By  evaluating  the  moment  of  vorticity  relative  to  a  line,  displaced  hy 
above  the  wall,  the  outer  flow  is  separated  from  the  inner  recirculation  region 
and  thus,  the  difficulty  of  determining  the  proper  ymax  Fmax  avoided. 

The  second  difference  between  (5.8)  and  (5.12)  is  the  substitution  of 
du/dy  for  the  |a;).  Since  the  derivatives  of  the  v  and  w  components  contribute 
very  little  to  the  vorticity  compared  to  the  du/dy  term,  the  vorticity  in  the 
boimdary-layer  model  is  approximated  by  du/dy.  By  making  these  changes, 
equation  (5.12)  provide  more  realistic  values  of  yma*  and  F^ax  for  the  separat¬ 
ed  flow  case. 

The  Klebanoff  intermittency  factor,  7,  given  by  (5.9)  must  also  be 
modified  to  take  into  account  the  displacement  of  the  shear  layer  from  the 
wall  by  the  distance  hy.  If  7  remains  the  s2Lme  function  of  y/S,  the  intermit¬ 
tency  is  too  high  at  the  outer  edge  of  the  shear  layer  in  the  separated  region. 
This  problem  can  be  avoided  by  letting  7=1.0  for  all  values  of  y  less  than  hy. 
hy  can  be  approximated  by  (yni«  —  b),  where  b  is  the  half  width  of  the  shear 
layer,  see  Figure  5.42,  and  is  estimated  for  separated  flows  from  the  following 
expression: 

U  _  ’T  ~  ^max)  ,  ,  .  fr  ■,n\ 

^  ~  A  p  (ymax  hy)  (5.13) 

’  C  max 


After  modifying  the  length  scale  of  (5.9)  to  zurcount  for  the  displzicement  of 
the  shear  layer  in  separated  flows  the  intermittency  factor  now  becomes: 

p  (y -Ymax + 

'“'kleb  — 


7  =  U.O  +  5.5 


2b 


(5.14) 


For  separated  flows,  the  intermittency,  as  defined  by  (5.14),  decreases  the 
eddy  viscosity  to  approximately  the  same  value  at  the  shear  layer  outer  edge 
as  it  would  have  achieve  at  the  outer  edge  of  an  attziched  boundziry- layer. 
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A  correlation  equation  for  the  experimental  profiles  in  the  separated  or 
recirculation  region  has  been  developed  in  order  to  estimate  the  moment  of 
vorticity.  This  procedure  was  necessary  because  of  the  limited  number  of 
measurements  within  the  shear  layer  which  makes  a  direct  determination  of 
du/dy  difficult.  The  correlation  equation  assumes  the  following  polynomials 
with  the  appropriate  boundary  conditions: 


u  =  ae(77)  +  (1  -  a)f(77)  +  cg(77) 

(5.15) 

where  u  and  77  are  given  by: 

U  -  Umin 

u  — 

(5.16) 

y  -  hy 

5-hy 

(5.17) 

Consider  the  following  polynomials  fimctions  e(r;),  f(77)  and  g(r/)  which  must 
satisfy  the  boundary  conditions  describing  the  separated  profiles. 


e(rf)  =  r/  -h 


(5.18) 


e^r/)  must  satisfy  the  boundary  conditions  at  ^=0,  e=0  and  at  77=1,  e=l  and 
e  =0.  Note  that  in  the  liinit,  as  rj  — »  0,  equation  (5.15)  is  linear  with  slope  of 
a.  The  f(  r? )  polynomial  is'  defined  as: 

f(fj)  =77^+277^-27?^  (5.19) 

with  the  boimdary  conditions  at  77=0,  f=0  and  at  77=1,  f=l  and  f'=0.  Thus, 
the  first  two  terms  on  the  right  hand  side  of  (5.15) ,  when  combined,  satisfy 
the  outer  boimdary  conditions  at  77=0.  The  final  polynomial,  g(77),  provides 
some  control  over  the  shape  of  the  profile  without  affecting  the  boundary  con¬ 
ditions  and  has  the  following  form: 

S(f7)  =r{^  -  277^  +  77**  (5.20) 


Note  that  g  must  satisfy  the  homogeneous  boundary  conditions,  at  77=0,  g=0, 
and  at  77=1,  g=0,  g'=0. 

The  polynomial-fitted,  separated  profiles  for  the  15  degree  step  zu'e 
presented  in  Figure  5.43.  The  agreement  between  the  curve-fits  and  the  ex¬ 
perimental  results  are  typical  of  the  separated  profiles  for  all  three  configura¬ 
tions,  so  only  the  15  degree  results  are  shown.  A  minor  problem  exists  with 
the  present  curve-fitting  method.  For  the  profiles  near  the  step,  the  slope,  a, 
becomes  very  large  which  causes  the  curve-fitted  polynomials  to  overshoot  the 
outer  boundary  conditions.  Thus,  the  boundary  conditions  are  satisfied  by  be¬ 
ing  approached  from  above.  This  is  undesirable,  but  only  complicated 
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transcendental  relations  are  able  to  avoid  this  oversnoot  problem.  The  curve- 
fitted  profiles  also  cam  not  approximate  the  reverse  flow  region  with  any  rea¬ 
sonable  degree  of  £u:curacy.  Although  there  are  some  problems  with  the 
present  curve  fitting  procedme,  the  resulting  polynomial  profiles  are  still  use¬ 
ful  for  the  evaluation  of  the  Baldwin-Lomax  parameters  because  it  is  only 
necessary  to  adequately  describe  the  middle  layer  of  the  profile,  not  the  region 
near  the  reverse  flow  or  at  the  outer  edge.  The  curve-fitted  profiles  shown  in 
Figure  5.43  illustrate  how  well  the  center  of  the  profiles  are  described. 

Comparing  and  C^y.  Ideally,  a  comparison  should  be  made 
between  the  shear  stress  computed  using  the  Baldwin-Lomax  method  and  an 
experimentally  measured  stress.  This  is  beyond  the  scope  of  the  present 
study.  Therefore,  an  indirect  procedure  has  been  used  which  is  based  on  a 
fundamental  assumption  in  formulating  the  Baldwin-Lomax  method.  It  is  as¬ 
sumed  that  the  Baldwin-Lomax  method  is  derived  from  the  empirical  observa¬ 
tion  of  Clauser  that  the  outer  layer  eddy  viscosity  is  determined  from: 

/it°  =  Kp4*UeTf  (5.21) 

where 

00 

4*  =  /  (I  -u/ue)dy  (5.22) 

0 

Since  the  vjelocity  profiles,  at  least  in  the  outer  layer,  have  been  meastired,  the 
value  of  4  cajfi  be  calculated  using  the  estimates  of  6  and  Ug  which  have  al¬ 
ready  been  discussed.  Equating  the  Baldwin-Lomax  formulas  for  /xt°  to 
Clauser’s,  (5.21),  the  constants  Cep  for  attached  flow  and  CepC^k  for  separat¬ 
ed  flow,  can  be  determined  from  the  following  relations: 


(5.23) 

~  F  v 

max  J  max 

4*U«Fn,ax 

(5.24) 

From  the  above  procedure  the  constants  for  the  Baldwin-Lomax  method  were 
determined  for  the  u  velocity  components  of  all  three  step  configurations  and 
are  presented  in  Figure  5.44.  The  dashed  lines  indicate  Baldwin-Lomax’s 
(1978)  recommended  results  for  the  two  constants.  The  comparison  of  Cep 
upstream  of  the  step  is  good,  with  an  average  result  of  approximately  1.4. 
Downstream  of  reattachment  the  average  value  of  C^p,  approximately  1.1,  is 
considerably  lower  than  the  1.6  suggested  by  Baldwin-Lomax  (1978)  for  at¬ 
tached  flow.  As  the  flow  continues  downstream,  C^p  continues  to  increase. 
This  result  suggests  that  at  10  step  heights  downstream  of  reattachment  the 
flowfield  has  not  totally  recovered  from  the  flow  separation  at  the  step. 
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Figure  5.44  Baldwin- Lomax  Constant  for  the 
Attached  and  Separated  Profiles 
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There  is  no  observable  effect  of  sweep  angle  on  Qp  upstream  of  the 
step.  The  results  seem  to  be  evenly  scattered  about  1.4  without  a  noticeable 
trend.  For  a  given  survey  location  downstream  of  reattachment,  the  effect  of 
increasing  the  sweep  angle  is  to  Increase  C^p.  These  results  indicate  that  as 
the  step  sweep  angle  increases,  the  flow  in  the  recovery  region  recovers  sooner 
from  the  effects  of  the  flow  separation  at  the  step.  This  is  in  agreement  with 
the  results  from  the  Coles’  law  analysis,  where  the  wake  pau’ameter,  IT,  de¬ 
creased  faster  with  increasing  step  sweep  angle  in  the  recovery  region. 

The  agreement  between  the  calculated  CepC^k  a^d  the  suggested 
results  of  0.4  is  very  good,  with  the  average  of  CepC^k  being  jiist  slightly 
higher.  The  0  degree  step  results  remainly  nearly  constant  at  a  value  of  ap¬ 
proximately  0.4.  The  results  for  the  15  degree  step  increased  slightly  just 
downstream  of  the  step  before  decreasing  at  a  x'/h  »2.5.  The  values  of 
CcpCvrk  ^he  30  degree  step  exhibit  a  trend  opposite  to  that  of  the  15  de¬ 
gree  step.  With  increasing  downstream  distance  from  the  step,  CepC^k  de¬ 
creases  until  a  location  of  x'/h  »2.5,  then  for  the  remainder  of  the  recircula¬ 
tion  region  the  values  increase. 

Correlation  of  Profiles  With  Baldwin-Lomax  Variables.  The  final 
analysis  was  to  examine  the  correlation  of  the  separated  profiles  of  both  the  u 
and  w  velocity  components  with  the  Baldwin-Lomax  variables.  Figures  5.45 
to  5.47  are  the  correlation  profiles  for  the  u  components  within  the  recircula¬ 
tion  region  for  each  step,  respectively.  The  profiles  are  plotted  with  respect  to 
a  nondimensional  velocity  scale  variable,  (u  —  Umijj)/(u  e  -  Umit.)  and  a  length 
variable,  (y  —  hy)/b.  The  0  degree  step  results,  Fi^e  5.45,  show  that  the 
sbe  separated  profiles  have  collapsed  into  a  single  curve  within  the  scatter  of 
the  data.  The  results  for  the  15  degree  step,  Figure  5.46,  show  slightly  less 
similarity.  This  is  because  the  results  near  the  step,  within  the  recirculation 
region  are  near  the  lower  limit  of  the  calibration  and  are  questionable.  The 
similarity  of  the  results,  shown  in  Figure  5.47  for  the  30  degree  step,  is  the 
poorest.  As  for  the  15  degree  case,  the  accuracy  of  the  30  degree  profiles  at 
x/h=0.3  to  1.0  and  below  y/h=1.0  are  questionable.  There  is  no  strong  trend 
evidenced  in  the  correlation  with  downstream  distance  using  the  Baldwin- 
Lomax  variables. 

The  correlations  of  the  w  velocity  components  using  the  Baldwin- 
Lomax  variables  are  shown  in  Figmes  5.48  and  5.49.  The  w  components  are 
all  equal  to  zero  for  the  0  degree  step,  so  only  the  results  from  the  15  and  30 
degree  steps  can  be  examined.  The  length  scale  used  for  the  w  component 
was  the  same  as  for  the  u  component  plots.  A  different  velocity  scale  was 
used  and  is  given  by  w/w^,  where  We=Ue  tan(A).  Figure  5.48  presents  the  w 
velocity  component  within  the  separation  region  of  the  15  degree  step.  The 
sunilarity  of  the  results  is  fair.  Figure  5.49  presents  the  correlation  of  the  w 
components  for  the  30  degree  step.  Again  the  similarity  is  fair.  One  notice¬ 
able  result  from  these  two  plots  is  that  the  ratio  of  w/wg  obtains  a  value 
Izirger  than  1.  For  the  15  degree  step  the  velocity  ratio  at  the  edge  of  the 
boundary- layer  is  approximately  1.25  and  for  the  30  degree  step  it  is  about 
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CORRELATION  OF  SEPARATED  PROFILES 
BALDWIN-LOMAX  VARIABLES 
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Figure  5.45  Correlation  of  the  u  Component  Separated 
Profiles  With  the  Baldwln-Lomax  Variable 
for  A  =  0  Degree 
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Figure  5.46  Correlation  of  the  u  Component  Separated 
Profiles  With  the  BaldMin-Loinax  Variable 
for  A  =  15  Degrees 
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Figure  5.47  Correlation  of  the  u  Component  Separated 
Profiles  With  the  Baldwin-Lomax  Variable 
for  A  =  30  Degrees 
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Figurie  5.48  Correlation  of  the  w  Component  Separated 
Profiles  With  the  Saldwln-Lomax  Variable 
for  A  »  15  Degrees 
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Figure  5.49  Correlation  of  the  w  Component  Separated 
Profiles  With  the  Baldwin- Lomax  Variable 
for  A  =  30  Degrees 
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1.15.  The  difference  seems  to  be  due  to  a  ttiming  of  the  free  stream  so  that 
the  velocity  is  directed  along  the  step.  This  is  contrary  to  the  flow  in  the  re¬ 
circulation  region  which  turns  so  that  the  w  component  is  flowing  in  the  nega¬ 
tive  direction,  along  the  step.  The  reason  for  this  is  not  clearly  understood. 

Coles’  law  and  the  modified  Baldwm-Lomax  variables  have  been  shown 
to  provide  techniques  for  correlating  the  data  from  this  study.  It  has  also 
be^  shown  that  the  flow  over  the  swept  backward-facing  step  wind  tunnel 
model  was  influenced  by  pressure  gradients  acting  in  the  flow  and  it  is  not 
clear  whether  these  gradients  are  due  to  initial  conditions  or  whether  the  step 
sweep  induces  them.  A  major  finding  from  the  Bzddwin-Lomax  analysis  was 
that  all  length  scales  used  in  the  separated  flow  region  must  take  into  accoimt 
the  displacement  of  the  shear  layer  from  the  wall  by  the  distance  hy . 
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Chapter  6 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


Summary 

The  objectives  of  studying  the  flowfield  behind  a  swept  backward¬ 
facing  were: 

(1)  to  provide  preliminary  mean  flow  measurements  behind  a  swept 
backward-facing  step. 

(2)  to  determine  the  magnitude  of  the  velocity  components  parallel 
and  normal  to  the  step  face. 

to  examine  the  effect  of  step  sweep  zmgle  on  the  velocity  field. 

(4)  to  examine  the  use  of  Coles’  law  of  the  wall  auid  wake  in  the 
attached  boundary-layer  region  of  the  present  swept  step  data. 

(5)  to  examine  the  use  of  a  modified  version  of  the  Baldwin-Lomax 
algebraic  turbulence  model  for  the  flow  over  a  backward-facing  step. 

(6)  to  correlate  the  spanwise  separated  flow  u  and  w  velocity 
components  using  the  Baldwin-Lomax  variables. 

Much  work  has  been  performed  on  the  flow  over  unswept  backward¬ 
facing  steps  but,  until  recently,  very  little  interest  has  been  shown  in  the 
swept  step  case.  Although  the  flo^^leld  is  quite  complex,  the  flow  over  a 
swept  backward-facing  step  is  one  of  the  simplest  three-dimensional 
reattaching  separated  flows.  This  type  of  flow  occurs  in  everyday 
engineering  problems,  such  as  the  flow  over  joints  in  swept  wings  of 
airplanes,  flow  around  buildings  and  over  clife,  and  the  flow  over  rotating 
bands  of  artillery  projectiles.  Also,  data  on  the  flowfield  over  a  swept 
backward-facing  step  may  be  useful  m  the  testing  of  turbulence  models  for 
romputational  fluid  dynamics  codes  because  it  is  a  unique  three-dimensional 
flow. 


The  present  study  was  performed  at  the  Chemical  Research 
Development  and  Engineering  Center,  Aberdeen  Proving  Ground.  The  flow- 
field  measurements  behind  a  0.5  inch  high,  backward-facing  step  with  sweep 
angles  of  0,  15  and  30  degrees,  were  conducted  in  the  Aerodynamics 
Research  and  Concepts  and  Assistance  Branch’s  subsonic  wind  tunnel.  A 
single  component  hot  film  anemometer  was  xised  to  measure  the  velocity 
parallel  and  normal  to  the  step  face.  The  component  of  velocity  normal  to 
the  surface  of  the  plate  was  assumed  small.  Oil  flow  observations  were  also 
made  in  order  to  determine  similarity  of  flowfields  with  the  studies  of  other 
investigators. 


Conclusions 


The  most  significant  conclusions  determined  from  this  study  are 
presented  below. 

(1)  The  oil  flow  observations  indicate  that  the  flowfields  over  the 
three  step  configurations  are  similar  to  those  obtained  by  other 
investigators.  Observed  reattaw:hment  distances  for  the  three  steps 
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axe:  0  degree  step,  x'/h  =  5.0;  15  degree  step,  x'/h  =  4.8;  auid  30 
degree  step,  x'/h  =  4.6.  From  the  oil  flow  tests,  the  reattachment 
length  decreases  with  increasing  sweep  angle. 

12)  The  mean  velocity  components,  u  and  w,  were  determined 
(assuming  v  to  be  negligible)  at  various  distances  up  and  downstream 
of  the  step  and  at  various  heights  above  the  step.  The  0  degree  step 
results  show  good  agreement  with  the  results  of  other  investigators. 

(3)  The  flow  in  the  recirculation  region  tiims  away  from  the  free 
stream  direction  until  it  is  approximately  normzd  to  the  step  face. 
Then,  at  reattachment  and  in  the  beginning  of  the  recovery  region 
the  flow  realigns  itself  with  the  free  stream  direction. 

(4)  From  the  hot  film  measurements,  as  the  step  sweep  angle 
increases,  the  flow  reattachment  distanced  decreases,  verifying  the  oil 
flow  results. 

(5)  The  boundary-layer  thickness  upstream  of  the  step  experiences 
typical  growth  for  a  turbulent  flowfield.  The  experimentally 
determined  6  at  the  step  is  in  excellent  agreement  with  the  result 
obtained  from  the  Schlichting  turbulent  boimdary-layer  thickness 
equation.  Within  the  recirculation  region,  £  decreases  and  remains 
nearly  constant  imtil  the  flow  enters  the  recovery  region.  Once 
downstream  of  reattachment,  S  increases  with  increasing  downstream 
distance.  Also  S  tends  to  increase  with  increasing  sweep  angle. 

(7)  Within  the  recirculation  region  the  y  component  stagnation  line 
is  approximated  by  a  second  order  polynomial.  As  the  step  sweep 
angle  increases,  the  magnitude  of  the  slope  of  the  line  also  increases. 

(8)  Downstream  of  the  step  the  w  component  stagnation  line  is 
approximated  by  a  straight  line.  These  resul*^s  indicate  that  at 
reattachment,  u  =  0,  the  w  component  continues  to  have  a  finite 
negative  value,  which  is  confirmed  by  the  oil  flow  observation.  As 
step  sweep  angle  increases  the  magnitude  of  the  slope  of  the  w 
component  stagnation  line  aJso  increases. 

(9)  The  u-component  velocity  profile  data  were  correlated  using 
Coles’  law  of  the  wall  and  law  of  the  wzdie.  A  least  squares  technique 
was  used  to- determine  the  four  constants,  Cf,  C,  H,  6  of  Coles’ 
profile  equation.  Upstrezim  of  the  step,  the  results  from  the  four 
constants  suggest  the  following:  for  the  0  degree  step  m  adverse 
pressure  gradient  is  present;  a  favorable  pressure  gradient  acts  on  the 
15  degree  step  flowfield;  for  the  30  degree  step,  a  slight  favorable 
pressure  gradient  acts  on  the  flow,  which  represents  a  nearly 
equilibrium  flow.  In  the  recovery  region,  the  skin  friction  coefficient, 
Cf,  generally  increases  with  increasing  downstream  distance  and 
sweep  angle.  The  intercept,  C,  increases  with  increasing  downstream 
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distance  iind  decreases  with  sweep  angle.  The  waJce  parameter,  IT, 
decreases  with  increasing  downstream  distance  from  reattachment, 
but  even  at  10  step  heights  the  value  of  IT  is  still  a  faotor  of  2  greater 
than  the  upstream  value.  Also,  11  decreases  with  increasing  sweep 
angle  indicating  that  the  larger  the  sweep  angle  the  faster  the  flow 
recovers  from  separation.  The  boundary-layer  thickness  determined 
from  the  Coles’  law  analysis  is  in  good  agreement  with  the  boundary- 
layer  calculated  directly  from  the  experimental  data  euid  summarized 
in  conclusion  (5) . 

(10)  The  semi-logarithmic  profiles,  at  the  step  for  the  three 
configurations,  show  evidence  of  the  adverse  pressiire  gradient  acting 
on  the  0  degree  step,  a  favorable  pressure  graidient  for  the  15  degree 
step  and  a  slight  favorable  gradient  acting  on  the  30  degree  flow.  The 
logarithmic  profiles  also  show  that  the  conditions  at  separation  for  the 
three  configurations  are  not  as  similar  as  the  y/h  versus  U/Ue  plots 
indicates.  The  logarithmic  plots  downstream  of  reattachment  for  the 
30  degree  step  shows  that  the  flow  is  recovering  from  the  separation 
with  increasing  downstream  distance. 

(11)  The  velocity  profiles  (u/ti^  versus  yjS)  determined  from  Coles’ 
law  ago  veil  with  the  experiment.  J  data.  The  profiles  upstream  of 
the  p*eo'  .ypically  collapse  into  a  single  ctirve.  Downstream  of 
reat.CL.v.ximent,  the  profiles  typically  become  fuller  as  the  downstream 

-stance  increases. 

(12)  The  applicability  of  the  Baldwin-Lomax  algebraic  turbulence 
model  to  the  results  of  this  study  was  also  investigated.  The  model 
was  modified  for  the  flow  over  a  bzu:kward-facing  step  by  accounting 
for  the  displacement  distance  of  the  she£ir  layer  from  the  wall  in  the 
separated  or  recirculation  region.  The  vzdue  of  C^p  upstream  of  the 
step  (approximately  1.4)  was  slightly  below  the  vadue  of  1.6  suggested 
by  Baldwin  and  Lomax.  Downstream  of  reattachment  the  value  was 
even  lower  (approximately  1.1)  but  increasing  with  increasing 
downstream  distance.  In  the  recovery  region.  Cep  tended  to  increase 
with  increasing  sweep  aingle.  Within  the  recirculation  region,  C<.p  C^t 
agreed  well  with  the  value  suggested  by  Baildwin  and  Lomax  of  0.40. 

(13)  The  separated  profiles  were  correlated  with  the  BaJdwin-Lomax 
variables.  Because  of  the  limited  number  of  data  points  describing 
the  separated  profiles,  the  profiles  were  curve  fitted  with  a 
polynomiaJ,  using  a  least  squares  technique.  This  was  necessary  in 
order  to  accurately  determine  a  value  for  du/dy,  which  was  required 
for  the  Baldwin-Lomax  method.  The  correlation  achieves  good 
results  when  the  displacement  of  the  shear  layer  from  the  wall  is 
taken  into  account.  The  u  component  profiles  indicates  good 
similarity  for  the  0  and  15  degree  steps  and  fair  similarity  for  the  30 
degree  step.  The  correlation  of  the  w  velocity  component  for  the  15 
and  30  degree  step  is  fair.  The  nondimensional  w  component,  w/w^. 
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achieves  values,  in  the  free  stream,  of  approximately  1.25  for  the  15 
degree  step  zind  1.15  for  the  30  degree  step. 

Recommendations 

Recommendations  for  future  resezu'ch  on  the  flow  over  swept 
backward-facing  steps  are  as  follows: 

(1)  More  detailed  measurements  downstrezim  of  the  step, 
especially  within  the  recirculation  region.  This  should  include 
more  survey  locations,  one  every  0.1  or  0.2  inches,  to  better 
define  the  reattachment  region  and  more  measurements  per 
profile.  The  measurements  for  each  profile  should  be  confined 
below  2.0  inches  (for  a  0.5  inch  step  height)  with 
measurements  every  0.010  inches. 

(2)  Experimentally  resolve  the  flow  reversal  ambiguity  in  the 
recirculation  region.  This  could  be  accomplished  by  several 
techniques  such  as  a  pulse  wire  anemometer,  split-film 
anemometer,  pressure  probe  or  shielding  and  rotating  a  single 
component  hot  wire  anemometer. 

(3)  Several  spauiwise  measurements  need  to  be  made  in  order 
to  confirm  the  SindependenceS  principle. 

(4)  Numerical  solutions  of  the  Navier-Stokes  equations  should 
be  comleted  for  the  flow  over  swept  backward-facing  steps. 
Further,  investigate  the  ability  of  computational  methods  to 
predict  the  flow  over  a  backward-feicing  step  and  the 
comparison  with  experimented  results. 
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LIST  OF  SYMBOLS 


a’s 
a,  c 


A+ 
A,  B 
AR 
b 
C 

Cep 


Cf 

Ckleb 

Cwk 


E 

F(y) 

Fjnax 

F(a) 

f{a) 

f(a  +  9) 
f(a  -  9) 


constants  used  in  reformulation  of  Coles’  Law 

constzints  used  in  polynomial  curve-fits  of 
separated  profiles  in  the  Bzildwin-Lomax  analysis 

Van  Driest  damping  function  =  26 

constzints  used  in  King’s  law 

aspect  ratio,  (width  of  wind  tunnel  model)  /h 

half  width  of  shear  layer 

intercept  of  the  semi-logarithmic  law  of  the  wall 

constant  in  Baldwin-Lomax  turbulence  model  for 
attached  flows  »=  1.6 

skin  friction  coefficient,  2(u,./u  g)* 

Klebanoff  intermittency  factor  =  0.3 

constant  in  Baldwin-Lomax  turbulence  model 
for  separated  flow  =  0.25 

output  voltage  from  hot  film  anemometer 

polynomial  function  used  for  polynomial  fit 
of  separated  profile  in  Baldwin-Lomax  analysis 

moment  of  vorticity 

maximum  of  moment  of  vorticity 

overall  velocity  ratio  based  on  a  -f- 

velocity  ratio  as  a  function  of  a 

velocity  ratio  as  a  function  of  a  +  5 

velocity  ratio  as  a  function  of  a  —  5 
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^iv) 

G{a) 

6{v) 

h 

hy 

K 

Lsep 

n 

Reji 

Re^ep 

Teal 

T 

Ta 

U 

Ud 

Ue 

Ueff 


polynomial  function  used  for  polynomial  fit 
of  separated  profiles  in  Baldwin-Lomax  analysis 

overall  velocity  ratio  based  on  a  —  6 

polynomial  fimction  used  for  polynomial  fit 
of  separated  profiles  in  Beddwin-Lomax  analysis 

height  of  the  step 

displacement  distance  of  the  shear  layer  from  the 
wail  in  separated  flows 

the  Clauser  constauit  =  0.0168 

distance  along  centerline  from  leading  edge  to 
top  of  step 

King’s  law  exponent  =  0.5 

Reynolds  number  based  on  step  height,  U  h/i/ 

Reynolds  number  based  on  separation  length, 

U  L^/u 

ambient  air  temperature  at  time  of  calibration, 

(•  F) 

ambient  air  temperature  at  time  of  velocity  profile 
survey,  ( ’  F) 

temperature  of  the  hot  film  sensor, 

(•  F) 

total  velocity 

V7u^  +  v^  +  w2)niax  -  V  +  w2)njj„ 

total  edge  velocity 

velocity  measured  by  hot  film  anemometer  at  some 
a  to  the  flow 

freestream  velocity 

component  of  velocity  perpendicular  to  step  face 
defined  by  u'*'  =  u/u 

velocity  scale  for  polynomial  fits  of  separated 
profiles  in  Baldwin-Lomax  method 
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Umax 

Umin 

Ut 

V 

w 

We 

X 

x' 

y 


ymax 


Z 

z' 

A 

n 

a 


7 

6 

V 


u-component  edge  velocity 

maximum  u  velocity  component 

minimum  u  velocity  component,  =0  for  attached 
flow,  negative  for  separated  flows 

shear  velocity,  y/r^/p  w 

component  of  velocity  perpendicular  to  wall 

component  of  velocity  paraJlel  to  step  face 

w  component  edge  or  average  velocity 

coordinate  in  streamwise  direction 

coordinate  normal  to  step 

coordinate  normal  to  wind  tunnel  test  section 
floor,  same  as  y' 

distance  above  surface  to  the  mziximum  of  the 
moment  of  vorticity 

defined  as  y  u^/t/ 

coordinate  normal  to  wind  tunnel  test  section 
floor,  same  as  y  coordinate 

lateral  coordinate  perpendicular  to  freestream 
direction 

coordinate  parallel  to  step  face 

step  sweep  angle 

Coles’  wake  parameter 

total  velocity  flow  direction  with  respect  to 
the  freestream  direction 

Klebanoff  intermittency  fewrtor 

boimdary  layer  thicknesses 

incompressible  displacement  thickness 

length  scale  for  polynomial  fit  of  separated  profiles 
in  Baldwin-Lomax  method 
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0 

K 

Mt 

A*t" 

V 

p 

T 

U) 


angular  offset  for  hot  film  measurements  =  15  degrees 

the  von  Karmaui  constant  =  0.43 

eddy  viscosity 

outer  region  eddy  viscosity 

kinematic  viscosity 

fluid  density 

shear  stress 

vorticity 


126 


